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METHOD AND APPARATUS FOR QUANTITATING SURFACE-BINDING 
OPTICAL RESONANCE PROFILES 
(HTS-016) 

5 FIELD OF THE INVENTION 

This invention relates to data quantitation methods and, in particular, to 
quantitation of data profiles obtained using optical resonance surface sensors. 

BACKGROUND 

1 0 Optical resonance surface sensors are utilized in a number of chemical, 

pharmacological, and biotechnology research fields, including proteomics and drug 
discovery. Quantitation of surface-binding optical resonance profiles, such as 
Surface Plasmon Resonance (SPR) profiles, allows real-time observation and 
analysis of molecular interactions, providing data uninfluenced by biochemical 

15 labeling. As a consequence, optical resonance surface sensors have become widely 
used for the characterization of biological surfaces and the real-time monitoring of 
binding events. 

At the present time, the most commonly used optical resonance surface 
sensors detect SPR, Surface plasmons are transverse electromagnetic charge-density 
20 waves that propagate parallel to the interface between a dielectric medium and a 
metallic film. Surface plasmons are generated by the interaction between the 
electron-rich surface of the metal and a charged particle or photon. Under 
appropriate conditions, the plasmons will resonate with light, resulting in the 
absorption of light. 

25 More specifically, at an interface between two transparent media of different 

refi-active index, light coming firom the side having higher refractive index is partly 
reflected and partly refracted. Above a certain angle of incidence, no light is 
refracted across the interface and total internal reflection (TIR) is observed. While 
incident light is totally reflected, the electromagnetic field component does penetrate 
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a short distance, on the order of tens or hundreds of nanometers, into a medium of 
lower refractive index, thus creating an exponentially detenuating evanescent wave. 
If the interface between the media is coated with a thin layer of metal, and the 
incident Hght is monochromatic and p-polarized, the intensity of the reflected hght is 
5 reduced at a specific incident angle. This produces surface plasmon resonance due 
to the resonance energy transfer between the evanescent wave and the surface 
plasmons. 

Many SPR sensors utilize the Kretschmann or Otto configurations, wherein 
the evanescent wave fi:om totally internally reflected monochromatic light traveling 

10 through a prism creates a surface plasmon in a metal film fliat is in contact with the 
material to be detected. In these configurations, the hght is incident on the metal 
fihn through the prism at an angle greater than a critical angle, known as the 
resonant angle ^spr- The surface plasmon is then detected by adjusting the angle of 
incidence until a sharp decrease in the reflected intensity is detected. If the 

15 permittivity or thickness of the dielectric layer changes, the resonant angle also 
changes. 

In another common SPR configuration, the light is incident firom the sample 
side of the metal film and is coupled to the surface to create the plasmon resonance 
by means of a diffiraction grating in the metal surface. In still another SPR 

20 configuration, the metal film is deposited on the outside of an optical fiber or other 
waveguide and hght is coupled into and out of the surface through the waveguide. 

The existence of the SPR is detected by measurement of the intensity of the 
reflected hght from the dielectric/metal interface, either from the back side 
(Kretschmann case) or the sample side (Otto and grating coupled cases). The 

25 resonance condition is sensitive to the effective refi-active index of the medium 
adjacent to the metal fihn, and hence to the configuration of its surface. The term 
effective refiractive index is used here because, if there is also a very thin (« lum) 
soUd biolayer in addition to the bulk fluid layer, then the biolayer will change the 
SPR angle and therefore cause the appearance of a change in the fluid refiractive 

30 index. This is the common case, where both the biofilm and the bulk fluid refiractive 
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index affect the signal. Essentially, everything in the evanescent layer a few 
hundred nanometers thick has an influence. The surface configuration is therefore 
also changed by any material adsorbed onto the metal film, so that the binding of 
biomolecules to the film resuhs in a chaQge of the effective refractive index of the 

5 dielectric. Becauseof this, adsorption ofmolecules on the metallic film or 

conformational changes in the adsorbed molecules can be accurately detected, SPR 
imaging can therefore be used to detect the presence and/or amount of a biopolymer 
on a chemically modified metal surface by quantitation of the change in the local 
index of refraction that occurs upon adsorption. 

10 For monochromatic or quasi-monochromatic illumination, the SPR angle will 

change directly according to the amount of bound materials; there is a linear 
relationship between the amount and the observed shift of the resonant angle. In 
particular, a linear relationship has been established between resonance angle and the 
mass concentration of biochemically relevant molecules such as proteins, sugars, and 

15 DNA. An SPR signal expressed in terms of angle shifts, or other units proportional 
to such shifts, is therefore a measure of mass concentration at the sensor chip 
surface. This means that analyte and ligand association and disassociation can be 
observed and rate and equilibrium constants may be calculated. 

In an alternative SPR configuration, a fixed angle of incidence is employed 

20 and the wavelength of the exciting light is varied, or a broadband optical source is 
enq)loyed together with an array spectrometer. With this approach, the SPR appears 
as a dip in the reflected intensity as a fimction of wavelength, with the wavelength of 
the SPR minimum varying linearly as a function of mass concentration at the chip 
surface. This configuration can be employed with any of the coupling mechanisms 

25 previously described (prism, grating, or waveguide). 

Surface plasmon resonance instruments may measure mass loading at a 
single spot, at several spots, or simultaneously at a large array of spots on an 
experimental stjrface (chip). In the case of an array instrument, the individual zones 
within which SPR shifts are measured are called Regions of Interest (ROIs). 



-3- 



PATENT 



The output of an SPR apparatus, whether the independent variable is angle or 
wavelength, is typically a graphed resonance curve. The location of the SPR 
resonance indicates the effective refractive index of the material on the sensor. 
Several approaches to SPR resonance curve shift quantitation have previously been 
5 utilized including, but not limited to: 

(1) 1st moment below a baseline 

(2) point of specific reflectance/signal 

(3) polynomial fit about the minimum 

(4) zero-crossing of the first derivative 
10 (5) reflectance/signal at a specific point 

(6) nonlinear fits of analytic functions 

Some literature distinguishes between Absolute Refi*active Index (RI) 
measurements and Relative RI measurements. The basis of this distinction is not 
entirely clear, since any precise absolute measurements must involve comparison 

15 with an RI standard. It is true, however, that, if a system is to be calibrated once and 
then used over time for absolute RI measurements, some algorithms are more 
suitable tiian others. This is due, for example, to such things as reduced sensitivity 
to sensor fouling or to baseline shifts. There also may be occasions where only 
changes in the refractive index are of interest. The same considerations apply, albeit 

20 to a reduced degree, to relative RI measurements of the sort required for measuring 
binding curves. 

The position of the SPR curve minimum is the most common indicator of the 
absolute refiractive index. The polynomial fit about the minimum (method (3)) and 
zero-crossing of the first derivative (method (4)) methods specifically determine this 
25 position. The other methods look at other aspects of the SPR curve in order to 
determine the minimum. Advantages of measuring absolute RI are that it looks at 
the SPR curve minimum, the measurement is insensitive to vertical shifts of the SPR 
curve (along the SPR signal axis), and the measurement is insensitive to fouling of 
the SPR sensor surface (which degrades the smoothness of the SPR curve)* 
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Disadvantages include that the measurement can be sensitive to the choice of points 
included in the analysis, potentially leading to distortion of the readmgs due to noise. 

Prior to performing any of the analysis methods, it is possible to smooth the 
SPR curve. This can lead to reduced noise when tracking refractive index versus 
5 time. One common smoothing algorithm is a least-squares smooth, which typically 
incorporates from 1 to 12 points on either side of a point when determinmg the 
smoothed value of that point Since smoothing does effectively bring multiple data 
points to bear on each smoothed point, algorithms which nominally depend on only 
one or a small number of data points may possibly be helped by smoothing 
10 procedures. In most cases, however, smoothing offers little or no benefit in well- 
designed fitting procedures, and can even hurt. 

(1) 1st Moment Below a Baseline Method 

The first moment method calculates the first moment of the SPR curve below 
a baseline. That is, only those parts of the SPR curve that are below the baseline are 
15 included in the calculation. Forthesimplecaseofnequally spaced data points, the 
algorithm may be expressed as: 

n 

- Baseline I *r 

1st moment = 

20 ~ " 

\SPRsignal , - Baseline \ 

f=i 

25 where the summations exclude all data points where the SPR signal is greater than 
the baseline. 

As long as the curve doesn't vertically shift, the algorithm will accurately 
track changes in the refractive index. Choosing where to set the baseline is not 
obvious. Frequently, the baseline is set at the midpoint of the SPR dip. Lowering 
30 the baseUne will include fewer points in the analysis and may lead to increased 
noise. Raising the baseline includes more points, but the resulting calculation will 
deviate even more from the curve minimum. If the points selected for the analysis 
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begin to shift out of the sensor range, the analysis will be degraded. There also may 
be times when an anomaly appears in the SPR curve, perhaps due to sensor surface 
fouling. It is still possible to use the sensor, so long as the anomaly is not near . 
the SPR dip. To deal with this situation, various regions may be selectively 

5 excluded from the analysis. Thus, even if the anomaly is below the baseUne, an 
accurate measurement results because the anomaly is excluded from the calculation. 

This technique can exhibit good performance with respect to shot noise or 
other similar random additive noises, in part because it involves a large number of 
points and a very simple algorithm. A helpful addition is the provision of 

10 interpolation calculations at the baseline cutoff levels on either wing to handle the 
fact that baselines will rarely pass directly through data points. Smoothing of these 
cutoff zones can also be helpful. However, the fundamental disadvantage of the 1^^ 
moment technique for SPR is its high sensitivity to intensity or signal basehne shifts 
(vertical shifts of the SPR curve) due to the inherent asymmetry of SPR resonances. 

15 High data point density is therefore required. This method also tends to be sensitive 
to sensor surface fouling. 

(2) Point of Specific Reflectance/Signal 

The point of specific reflectance/signal method utilizes the pixel position at 
which the SPR curve is a pre-specified value. The curve is initially examined in 

20 order to approximately locate the data point where the curve is closest to the 

specified value. An nth order polynomial least squares fit is then performed in order 
to interpolate and identify exactly where on the curve the value occurs. The search 
may be done from either the left or the right side of the curve. As the SPR curve 
shifts along the x-axis, the identified point will follow this shift. 

25 This method is susceptible to alterations in the overall shape of the SPR 

curve and it does not directly determine the minimum point. Still, it is useful for 
quantitating small shifts. In addition, it can be utilized to expand the dynamic range 
of the sensor. For example, if the minimum of the SPR curve is below pixel #1 (off 
the left side of the sensor range) the sensor may still be employed to track some 

30 other point on the SPR curve. 
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This technique has little to recommend it except simplicity. It suffers from 
much greater sensitivity to intensity or baseline changes than method (1), as well as 
being very sensitive to small changes in resonance shq)e. It uses very little of the 
resonance curve data for each determination, and thus has poor noise transfer 
5 performance. For decent performance, data point spacing must be tight. 

Another version of this technique is to perform the calculation on both sides 
of the resonance and use the mean of the two to track the resonance. This is then 
similar to method (1), but with poorer noise performance since, again, only a 
relatively small subset of the data is used. It also has greater sensitivity than method 
10 (1) to intensity shifts. 

(3) Polynomial Fit About the Minimum 

In the polynomial fit about the minimum method, a first pass of the SPR 
curve is made in order to find the approximate location of the minimum of the 
resonance. An nth order polynomial least squares fit (using some number of points 

1 5 about the minimum) is then performed in order to interpolate the position of the 
minimmn. This method is suitable for absolute refiractive index measurements. This 
method is not susceptible to y-axis shifting of the SPR curve. However, siace it only 
uses a relatively small number of points on the curve, it can be susceptible to the 
selection of points to be included in the calculation. This can sometimes lead to 

20 occasional anomalies m the analysis results (noise). 

This method works fairly well, especially with respect to intensity and 
baseline shifts, but depends on relatively little data and so exhibits sub-optimal noise 
transfer performance. As with any method that actually attempts to locate the 
minimum of the resonance (as opposed to its overall position), it can be sensitive to 

25 noise near the bottom. Higher order polynomials, which can better fit the shape of 
the resonance, allow a wider range of data points to be included, but also offer more 
opportunity to move the minimum around to accommodate noise on individual data 
points. In other words, higher order polynomials have excess degrees of freedom 
that translate into increased noise injection into the shift determination. In practice, 
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tight data point spacing is required, but the needed angle scan range can be fairly 
limited, 

(4) Zero-Crossing of the First Derivative 

The zero-crossing of the first derivative method is derived from the fact that 

5 the first derivative changes sign about the minimum of the SPR curve. In the zero- 
crossing algorithm, the approximate minimum point of the SPR curve is initially 
determined, a linear least-squares fit of the fu-st derivative is performed using a few 
points about the approximate minimum, and the zero-crossing point is interpolated. 
This amounts to a more or less direct determination of the actual minimum of 

10 the resonance curve, which is perhaps the least well-defined point on the curve. 
Many algorithms can be used to estimate the derivative and to find its zero crossing, 
and most of these do include some inherent smoothing or multi-point fitting, so that 
mathematically this process can be equivalent, for example, to method (3) 
(polynomial fit about the minimum). Using Savitsky-Golay derivatives, for 

15 example, the process is fixlly equivalent to method (3). Depending on the detailed 
implementation, this method can perform as well as method (3), and has much the 
same advantages and disadvantages. Poorly implemented, it can exhibit very poor 
noise performance. Again, tight data point spacing is needed. This method is 
suitable for making absolute refractive index measurements. 

20 (5) Reflectance at a Specific Point 

The reflectance/signal at a specific point method does not follow shifts of the 
SPR curve. It merely looks at the value of the SPR curve at a particular angle (or 
wavelength) position. If necessary, an nth order polynomial least squares fit is 
performed in order to interpolate the signal value at the desired Specific Point. This 

25 technique has often been used in the academic literature, since it can be utilized in a 
system without either moving parts or an array detector. It depends on the linearity 
of the side-wall of a resonance over a modest range of SPR shifts. It is essentially 
equivalent to method (2) (point of specific reflectance/signal) for small shifts, but, 
unlike method (2), it cannot accommodate larger shifts. One described 

30 implementation uses array detector information in which several data points 
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surrounding the Specific Point are fitted, thereby somewhat improving the noise 
performance. 

(6) Nonlinear Fits of Analytic Functions 

Theoretical Responses. It has been common in the academic hterature of 
5 SPR sensors to compare measured angular response curves with theoretical 

responses, and in some cases fits to these theoretical curves have been performed as 
a means of measuring fihn thicknesses, which amounts to measuring resonance 
shifts. 

A paper from Wolfgang Knoll's laboratory in Mainz [M. Zizlsperger and 

10 Wolfgang Knoll, "Multispot parallel on-line monitoring of interfacial binding 
reactions by surface Plasmon microscopy", Progr. Colloid Polym. Sci. 109: 244-253 
(1998)] shows diagrams (Figures 2 and 5) which appear to represent fits of 
experimental SPR angle response curves to simple Fresnel multilayer theory. The 
fits are fahrly poor, but would suffice to measure resonance shifts. There is no 

15 discussion of the fitting method. Possible methods include manual trial and error 
parameter adjustment, or non-linear methods such as Levenberg-Marquardt. 

A paper on Grating SPR sensors [C,R. Lawrence, N J, Geddes, D,N. Furlong, 
J, R. Sambles, "Surface Plasmon resonance studies of immimoreactions utilizing 
disposable diffraction gratings", Biosensors & Bioelectronics 11: 389-400 (1998)] 

20 shows fits of experimental data to GSPR theory using a coupled wave method [J. 
Chandezon, M.T. Dupuis, G. Comet, and D. Maystre, "Multicoated gratings: a 
differential formalism applicable in the entire optical region", JOSA 72: 839-846 
(1982)]. The fitting procedure is not described, but could be a conventional non- 
linear least squares fit. It is stated that the fit involves a large number of parameters, 

25 including the optical constants of gold, the pitch and groove depth of the grating, and 
a groove distortion factor, as well as the thickness and (sometimes) the optical 
properties of an adsorbed protein layer. In some cases, certain of these parameters 
may be fixed. The fits are very good, indicating that great care was taken 
experimentally to avoid distortions and artifacts. 
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These techniques are in principle capable of very good results for SPR shift 
quantitation, provided that the number of free parameters is kept small. In other 
words, if a preliminary fit were used to establish all of the parameters except for the 
adlayer thickness, and subsequent run-time fits left only the thickness and a couple 
5 of amplitude scaling parameters free, then performance should be reasonably good. 
Some downsides to this approach are: 

a) The model calculations are very complex and intensive, particularly 
fertile GSPR case. 

b) Nonlinear fits tiiemselves are very compute intensive (beyond the 
10 necessary theoretical model evaluations embedded therein), and exhibit instabilities 

that complicate automated operation. 

c) In practice, frequent deviations of resonance shapes from the 
theoretical ideal are observed. Indeed, the data of Zizlsperger et al show rather poor 
fits. 

1 5 Explicit Functions. Yet another approach that has been suggested is non- 

linear fitting of relatively tractable explicitly defmed mathematical functions to the 
run-time data. Unlike the Theoretical Responses, which are numerical profiles 
resulting from a complex computational process for a given set of input parameters, 
such explicit functions could be computed cheaply and are therefore much more 

20 attractive. In a way, this is an extension of flie polynomial fitting approach, although 
polynomials have the property of allowing fast linear least squares fits. 

If functional forms could be identified that have the flexibility to fit well to 
the range of observed profile shapes, this approach can work well. The downsides 
are: 

25 a) The difficulty of identifying such functional forms, unless we restrict 

the fits to a small region around the minimum, in which case see Polynomial Fits, 
above. As a corollary, the difficulties of ensuring that future oddities in profile 
shapes fliat are not modelable by the chosen functional form do not occur. 

b) The need, in general, for non-linear least squares or otiier non-linear 

30 fitting procedures. 
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c) Degraded noise performance due to the extra shape parameters Ukely 
needed to represent the full range of observed profiles. This problem can be 
sidestepped by the approach suggested above, in which a preliminary fit is used to 
fix most of the fi:ee parameters, leaving only resonaace shift and amphtude(s) as 

5 fitting parameters. 

Although most of the foregoing discussion has been about the most 
commonly utilized form of surface-binding optical resonance, SPR, optical 
resonance curves obtained fi'om other types of surface-binding optical resonances are 
quantitated in a similar fashion. Examples of devices and techniques for measuring 

10 such resonances include optical waveguide sensors (such as the BIOS-1 angle 
scanned grating coupled optical waveguide instrument currently made by Artificial 
Sensing Instruments AG), grating couplers (K. Tiefenthaler, "Grating couplers as 
label-free biochemical waveguide sensors'*, Biosensors Bioelectron., 8 (1993) xxxv- 
xxxvii), plasmon waveguide resonance devices (Z. Salamon, H. A. Macleod and G. 

1 5 Tollm (1997) Coupled Plasmon-Waveguide Resonators: A New Spectroscopic Tool 
of Probing Film Structure and Properties. BiophysicalJoumal 73:2791-2797), 
diffraction anomaly sensors (U.S. Pat. No. 5,925,878, Challener, 2000, "Diffraction 
anomaly sensor having grating coated with protective dielectric layer"), resonant 
mirror devices (R. Cush, J. M. Cronin, W. J. Stewart, C. H. Maute, J. MoUoy and N. 

20 J. Goddard, "The Resonant Mirror: a Novel Optical Biosensor for Direct Sensing of 
Biomolecular Interactions", Biosensors & Bioelectronics 8(7/8), 347-353 (1993)), 
long range SPR (http://plazmon.ure.cas.cz/tobiska/optsen01.pdf), and the Perkin 
Elmer optical resonance analysis system described in WO9930135 (Tracy et al, 
published June 17, 1999), To some extent, these devices and techniques are 

25 variations on two themes - dielectric planar waveguides and metal film SPR, or 
combinations thereof. All require no fluorescent label and are often called 
biosensors, but they are not restricted in their ^pUcation to characterization of 
biofilms. 

What has been needed, therefore, is a way to quantitate surface-binding 
30 optical resonance curves that requires low computational complexity and a minimal 
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number of scan data points while still providing accurate determination of the 
resonance angle and good noise performance. This ability will in turn allow for 
longer scan times and/or scans over a greater number of ROIs, enhancing the utility 
of optical resonance surface sensor techniques for the observation of, for example, 
5 such phenomena as adsorption onto chemically modified metal surfaces, binding 
events involving biological molecules such as DNA, proteins, enzymes, and 
antibodies, and immunologic phenomena such as antigen-antibody reactions and 
antigen stimulation of tissue. 

10 OBJECTS OF THE INVENTION 

Accordingly, an object of the present mvention is to provide a new 
mechanism for accurate quantitation of surface-binding optical resonance curves. 

A particular object of the invention is to provide a method for quantitation of 
surface-binding optical resonance curves that has relatively low computational 
15 complexity. 

A further particular object of the present invention is to provide a method for 
quantitation of surface-binding optical resonance curves that requires a minimal 
number of scan data points. 

Anodier particular object of the invention is to provide a mechanism for 
20 quantitation of surface-binding optical resonance curves that is less sensitive to shot 
noise. 

Yet another particular object of the invention is to accommodate the use of 
long surface-binding optical resonance scan times. 

Still another particular object of the invention is to provide a mechanism for 
25 quantitation of surface-binding optical resonance curves that provides increased scan 
speeds. 

Another particular object of the invention is to provide the ability to obtain 
surface-binding optical resonance results for a sample more quickly than is presently 
possible. 

30 Still a further particular object of the invention is to accommodate surface- 
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binding optical resonance scans over a greater number of ROIs. 

A further particular object of the mvention is to provide vaUd extension of 
the dynamic range of surface-binding optical resonance scans, so that useful results 
may be obtained even when only partial resonance curves are available due to 
5 hmitations on the angular scanning range employed. 

SUMMARY 

These and other objectives are met by the present invention in which 
Empirical Profile Fits (EPF) are used to quantitate surface-binding optical resonance 

10 profiles. The EPF method and apparatus use observed surface-binding optical 
resonance profiles at a relatively fine angle or wavelength spacing to create an 
empirical functional form. In the preferred embodiment, surface-binding optical 
resonances to be quantitated are then fit by a three parameter model consisting of (a) 
the empirical profile itself, (b) a numerical derivative of the profile generated using 

15 Savitsky-Golay filters, and (c) a constant additive offset to model any baseline shifts 
or continuum illumination. The shift of the fitted resonance is determined as the 
ratio of the derivative fit coefficient to the profile fit coefficient. When the shift so 
determined exceeds a predetermined threshold, the original EPF profile is shifted by 
Lagrange Interpolation to the approximate location of the resonance, and the fit is re- 

20 done. 

The EPF peak finding process has two stages, the Calibration stage and the 
Fit stage. In the Calibration stage, a calibration scan is obtained with relatively fine 
angle or wavelength spacing over a range sufficient to include full resonance profiles 
for all regions. Smoothed, subsampled model profiles for each Region of Interest 

25 (ROI) are generated, together with first derivative curves and diagnostic information. 
In the Fit stage, individual ROI scans are used for measurement of resonance shifts 
relative to the calibration model. The run-time scans may be at much coarser angle 
or wavelength spacing, may cover entirely different angle or wavelength ranges than 
the calibration data, may have either scan polarity and need not encompass the mtire 

30 surface-bmding optical resonance profile. Estimates of peak position can be made 
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even from fragmentary data on one side of the resonance, so that useful resuUs may 
be obtained even when only partial curves are available. A minimum of three data 
points is needed for each determination, although normally more points are used. 
Resonance shifts, estimated absolute angles or wavelengths, time of resonance 
5 minimum, and additional diagnostic and quality information are retumed. Data 
obtained from either calibration or fit runs may optionally be exported for analysis 
using other systems. 

In the preferred embodiment, the two stages of the process involved in EPF 
peak finding are implemented via two separate modules in software, the Calibration 

10 module and the Fit module. When a new chip is put into service, or a new 
calibration is desired on an old chip, a calibration scan is performed. During the 
calibration scan, surface-binding optical resonance data are obtained under static 
conditions over a relatively long angle or wavelength range, v^th relatively small 
angle or wavelength steps. The raw scans obtained during the calibration scan, one 

15 for each ROI being calibrated, are then optionally smoothed and subsampled to 
increase the point density. In the preferred embodiment, Lagrange or spline 
interpolation is used to inq)lement the subsampling. Derivatives of the subsampled 
profiles are computed, using a Savitsky-Golay filter or any other suitable method 
known in the art. Properties of the profiles are then calculated, mcluding, but not 

20 limited to, approximate resonance position, depth, and width. These results are then 
stored in preparation for later fitting procedures. 

In the Fit stage, angle or wavelength scans are measured for particular ROIs 
during an experiment. These scans may be obtained at much lower point density 
than the calibration scans, possibly over shorter angle or wavelength ranges. The role 

25 of the fitting module is to look at these experimental scans, identify the region 

encompassing the resonance, and fit that region using the previously stored empirical 
profiles in order to quantify and return the desired values, including the shift in the 
resonance as compared to its angular location at calibration time. While absolute 
angular positions or wavelengths generally cannot be measured accurately, since the 

30 resonance position is itself not well defined, rough values can be obtained by adding 
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the shifts to the nominal peak location determined during calibration, and these are 
also returned. 

In a preferred embodiment of the EPF calibration and fitting system, the 
Instrument Control and Data Acquisition Software gets and sets internal parameters 
5 in the EPF Calibration module and then sends raw data from a calibration scan to the 
EPF Calibration module, which funnels the data through a Subsampler and a 
Savitsky-Golay smoothing routine before taking derivatives and characterizing the 
data to create a model calibration profile for the chip being scamied. The model 
calibration profiles are then optionally stored. Next, the Instrument Control and 

10 Data Acquisition Software gets and sets internal parameters in the EPF Fitting 
module and then sends raw data from a run-time scan performed utilizing the chip to 
the EPF Fitting module, which qualifies it, queries the EPF Calibration module for 
the model profile for the chip, and then fits the curve utilizing various matrix 
routines, iterating the process when necessary. Results from the fitting process are 

15 then returned to the Instrument Control and Data Acquisition Software, which 
provides them to the user. 



BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 illustrates the steps of an embodiment of the method for quantitating 
20 surface-binding optical resonance profiles according to the present invention; 

Fig. 2 is a block diagram of an example implementation of the apparatus for 
quantitating surface-binding optical resonance profiles of the present invention; 

Fig. 3 is a diagram illustrating the graphical output of a calibration scan 
according to the present invention; 
25 Fig. 4 is a diagram illustrating the graphical output of a run-time scan 

according to the present invention; 

Fig. 5 is an operational flowchart of the calibration stage of the present 
invention; 

Fig. 6 is an operational flowchart of the processing of a calibration profile; 
30 Fig. 7 is a diagram depicting the definitions of the profile properties returned 
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during the calibration stage of the present invention; 

Fig. 8 is a screenshot depicting an example user interface for the initiation of 
the calibration stage; 

Fig. 9 is a screenshot depicting an example output from the calibration stage; 
5 Fig. 1 0 is a screenshot depicting a different example output from the 

calibration stage, wherein there is a "bad ROI" 41 showing no resonance; 

Fig. 1 1 is a screenshot depicting yet another example output from the 
caUbration stage, showing the details of "good ROI" 190; 

Fig. 12 is an operation flowchart for the optional chip qualification 
10 procedure; 

Fig, 13 is an operational flowchart of the fitting stage of the present 
invention; 

Fig. 14 is an operational flowchart of the processing of one ROI scan; 
Fig. 15 is an operational flowchart of the determination of the resonance 
1 5 angle utilizing an empirical profile fit; 

Fig. 16 is a screenshot depicting an example output from the fit stage; 
Fig. 17 is a screenshot depicting an example embodiment of a local error log, 
showing reports of bad fits to particular ROIs; and 

Fig. 18 is a screenshot depicting an example embodiment of an error log on a 
20 remote networked supervisory computer, showiag reports of bad fits to particular 
ROIs. 

DETAILED DESCRIPTION 
The current invention utilizes Empirical Profile Fits (EPF) to quantitate 
25 surface-binding optical resonance profiles. The EPF peak finding process has two 
stages, the Calibration stage and the Fit stage. In the CaHbration stage, a calibration 
scan is obtained with relatively fine angle or wavelength spacing over a range 
sufficient to include fiiU resonance profiles for all regions. Smoothed, subsampled 
model profiles for each Region of Interest (ROI) are generated, together with first 
30 derivative curves and certain diagnostic information. In the Fit stage, individual ROI 
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scans are used for measurement of angle or wavelength shifts relative to the 
calibration model. The run time scans may be at much coarser angle spacing, may 
cover entirely different angle ranges than the calibration data, and may have either 
scan polarity. Resonance shifts, estimated absolute angles or wavelengths, time of 

5 resonance minimum, and additional diagnostic and quality information are returned. 
Some precursors of certain components of the empirical profile fit of the 
present invention are disclosed in two prior patents. The use of derivatives in a fit 
to more accurately calculate and correct for wavelength shifts is disclosed in U.S. 
Pat No. 5,303,165 (Ganz et al, 1994) at columnl 1, line 1 ; column 16, Une 40 and 

10 U.S. Pat No. 5,308,982 (Ivaldi et al, a994) at column3, line 27 thru columnS, Une 
43 . In addition, the general idea of using empirical profiles for fitting has been 
disclosed for use in other contexts. For example, while the words "empirical profile" 
are not used, a somewhat similar idea is described in U.S. Pat. No. 5,308,982 at 
column 3, line 46, as a model matrix of known form. 

1 5 Empirical Profile Fits have been experimentally shown to work at least as 

well as or better than, polynomial fits in quantitating Surface Plasmon Resonance 
profiles. Shot noise limited performance is typically a few tens of microdegrees. 
EPFs have the further advantage of allowing much sparser angular point spacing, 
which can increase scan speed several-fold. This is especially true when long scans 

20 are needed to acconmiodate a wide range of resonance shifts in certain experiments. 
EPF SPR quantitation can be used with any implementation of a mechanism for 
creating and measuring Surface Plasmon Resonances, including, but not limited to, 
diffraction grating, prism, and waveguide implementations, with either angle or 
wavelength as the scanning variable. It may also advantageously be used to 

25 quantitate shifts of other resonances used for chemical or biochemical surface 
sensors, such as waveguide sensors and resonant cavity sensors. 

For simplicity, the description of the preferred embodiment that follows 
specifically discusses the parameters used when quantitating results obtained 
utilizing a grating-coupled angle scan SPR device. However, the EPF quantitation 

30 technique of the present invention is not limited to quantitation of data obtained 
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utilizing SPR devices, nor is it limited to quantitation of data obtained utilizing angle 
scan techniques. The EPF quantitation technique of the present invention may be 
advantageously employed to quantitate data produced using any of the many 
varieties of surface-binding optical evanescent wave sensor technologies which 
5 employ angle scanning, wavelength scanning, or a combination thereof (such as 
predispersion or correlated scanning) and which produce resonances similar to SPR 
resonances that shift in angle or wavelength more or less linearly with surface 
adlayer mass or binding. Examples of techniques producing resonances suitable for 
quantitating using EPF include, but are not limited to, the various techniques 

10 previously discussed in the Background section. 

Some of these techniques may produce sharper or broader resonances than 
SPR, the resonance position may be more or less sensitive to biolayer mass density, 
and/or they may exhibit different resonance propagation lengths, allowing for higher 
or lower ROI densities in array configurations. The EPF procedure still may be 

15 applied directly to quantitate resonances produced using any of these techniques. 
The main difference in application from the preferred embodiment described herein 
is that the size of angle steps is scaled to the width of the resonances (or, in the less 
frequently used wavelength mode, the FWHM of the resonances in wavelength 
units). In addition, while the curves obtained using some of these techniques exhibit 

20 positive peaks rather than negative-gohig dips as resonances, this has no effect on 
EPF except with respect to the profile acceptance criteria. Further, while various of 
these sometimes use a label and use the evanescent wave to excite fluorescence, so 
long as the exciting beam is angle scanned, EPF can still be used, even though 
diffiise fluorescence is being measured rather than reflected intensity. 

25 In general, the window over which the fit is performed can be adjusted at 

will, but is usually chosen to cover the main part of the resonance dip, leaving out 
the far wings. For best results, resonances at each ROI are fitted using Empirical 
Profiles measured for that ROI. However, good results are also obtained using 
Empirical Profiles from other ROIs exhibiting similar resonance shapes. The noise 

30 performance of EPF is generally very similar to that for analytical profile fits 
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employing the same more or less orthogonal parameters: intensity, resonance shift, 
and baseline offset. 

The EPF procedure does not attempt to measure the absolute location of a 
resonance, nor does it even include any definition of "resonance center" . Rather, it 

5 enables precise measurement of shifts or changes of resonance position from the 
position of the original model resonance. However, in order to simplify 
interpretation, an "absolute angle" is arbitrarily assigned to the initial model profile, 
simply by finding a measure of its minimum as the minimum of a parabolic fit to the 
bottom n points. Typically, n=5. This semi-arbitrary assignment allows reporting of 

10 "absolute angles" which are in the ballpark of other measures. The data points used 
in the EPF fit need not be at the relatively fine angular spacing used during the 
original calibration scan. In general, any data point spacing may be used which 
allows at least three points to be found somewhere within the resonance. Using 
coarser spacing may, however, have an impact on the precision of determination of 

15 the resonance angle. 

Overview, The currently preferred embodiment of a system for 
implementation of the present invention is an array system in which a 2-D CCD 
detector is used to image a sample chip containing 100-400 or more ROIs. Quasi- 
monochromatic light fi:om a filtered, p-polarized LED (Light-Emitting Diode) is 

20 used for illumination. The angle ofincidenceofthe light beam on the chip is varied 
by mechanically scanning the light source. For each angle of incidence, one or more 
CCD frame exposures are used to perform a simultaneous measurement of the 
reflected light intensities for all the members of the ROI array. While the preferred 
embodiment of the EPF system of the present invention is specifically targeted 

25 toward such array systems, it works equally well on single or few chaimel 
instruments. 

In this embodiment, when a new chip is put into service, or at any other time 
desired, a calibration scan is obtained with relatively fine angle spacing over an 
angle range sufficient to include fixll resonance profiles for all regions. The 
30 preferred angle range will depend on the refractive index (RI) of the buffer, any solid 
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films or adlayers present at the sensor surface, the sensor material, and (if 
applicable) the grating period. For example, when using grating coupling, a gold 
metal film, aqueous samples, a center wavelength of 875nm filtered by an 
Interference Filter to ~4nm F WHM, and an 840 nm grating groove spacing, the 

5 preferred angle range is typically 20° ± 2 In a preferred embodiment of the 
invention, the angle is controllable to millidegrees and can be measured to 
microdegrees. All these parameters can be adjusted, resulting in different nominal 
angles and possibly different FWHM of resonances, creating a need for different 
step sizes. In addition, if Wavelength scan mode is used, step sizes are in 

10 wavelength. To convert the exan:q)le above to Wavelength scan mode directly, the 
angle is fixed at -'20 deg and wavelength is scanned roughly from 860 to 890nm. 
The resonances would be - 10 nm wide. Similarly, for other Optical Resonance 
systems, the nominal wavelength, angle, and FWHM might be very different, 
requiring scaling of the various EPF parameters, but the process remains identical. 

15 In flie preferred embodiment, the calibration data set is submitted to the main 

calibration procedure, epfCal, which generates smoothed, subsampled model profiles 
for each ROI, together with first derivative curves and certain diagnostic 
information. This calibration set is kept in memory, and may optionally be archived 
for future use. Angle step sizes need not be uniform, as actual measured angles are 

20 supplied to the main calibration procedure for each data point. 

At run time, as experimental kinetic or other data are being obtained, 
individual ROI scans are submitted to the main run time fit procedure, epflFit, for 
measurement of angle shifts relative to the calibration model. The run time scans 
may be at much coarser angle spacing, and may cover entirely different angle ranges 

25 than the calibration data. They may also have either scan polarity (i.e., either up or 
down angle). Actual measured encoder angles must be supplied for each data point, 
and times of measurement of each data point must also be provided. Angle shifts, 
estimated absolute angles, time of resonance minimum, and additional diagnostic 
and quaUty information are returned. 

30 The run time scans need not encompass the entire SPR profile. Estimates of 
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peak position can be made even from fragmentary data on one side of the resonance, 
A minimum of three data points is needed for each determination, although normally 
more points are used. This ability to provide valid "edge info" for SPR scans, so that 
useful SPR results may be obtained even when only partial curves are available, is 

5 one of the many advantages of the present invention. A further advantage may be 
had by the optional ability to export data obtained from either calibration or fit runs 
for re-analysis, or additional analysis, using one or more other analysis systems. 

The process involved in EPF peak finding naturally divides into two parts, 
which, in the preferred embodiment, are implemented via two separate modules in 

10 software. These are the Calibration stage and the Fit stage. As seen in Figure 1, 
when a new chip is put into service 1 10, or a new calibration is desired on an old 
chip, a calibration scan 120 is performed. During the calibration scan 120, SPR data 
are obtained under static conditions over a relatively long angle range, with 
relatively small angle steps. The choice of what fluid is present on the sample side 

15 of the sensor surface depends on what measurements are being performed. The 
angle range is then chosen to accommodate the refractive index of the fluid. The 
precise values of these parameters need not be fixed, and can be varied as experience 
demands, but nominal 0.05 deg steps, as employed in the preferred embodiment, are 
highly suitable. These steps are refened to as Base Steps and are indexed by Base 

20 Index values. Actual angle values obtained from the angular encoder accompany 
each intensity value. Units for system angles may be degrees or other units. Any 
linear units are suitable, so long as they are used consistently. 

The raw scans obtained during the calibration scan 120, one for each ROI 
being caUbrated, are then optionally smoothed and subsampled 130 to increase the 

25 point density by N, where N typically equals 1 to 10. In the preferred embodiment, 
Lagrange or spline interpolation is used to implement the subsamplmg, although any 
other method known in the art would be suitable. Derivatives of the subsampled 
profiles are computed, using a Savitsky-Golay filter or any other suitable method 
known in the art. In the preferred embodiment, certain properties of the profiles are 

30 then optionally calculated, including, but not limited to, approximate resonance 
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position, depth, and width, in order to qualify the resonances as adequate and 
correct. These results are then stored 140 in preparation for later fitting procedures. 

In the Fit stage, angle scans are measured for particular ROIs during an 
experiment 150. These scans may be obtained at much lower point density than that 

5 of the calibration scans, and possibly over shorter angle ranges. In the preferred 
embodiment, the actual encoder angle and actual midpoint-of-integration time for 
each point are provided to the fitting module 160. The role of the fitting module 160 
is to look at the experimental scans obtained during the experiment 150, identify the 
region encompassing the resonance, and fit that region using the previously stored 

10 empirical profiles 140 in order to quantify and return 170 the desired values, 

including the shift in the resonance as compared to its angular location at calibration 
time. While absolute angular positions generally cannot be measured accurately, 
since the resonance position is itself not well defined, rough values can be obtained 
by adding the shifts to the nominal peak location determined during calibration, and 

15 these are also returned 170. 

Fig. 2 depicts a block diagram of a preferred embodiment of the EPF 
calibration and fitting system, including links to the instrument control software. As 
seen in Fig. 2, the Instrument Control and Data Acquisition Software 202 is utilized 
to get and set internal parameters 204 in the EPF Calibration module 210. In this 

20 embodiment, raw data 212 firom a calibration scan is sent fix)m the Instrument 
Control and Data Acquisition Software 202 to the EPF CaUbration module 210, 
which fimnels the data through a Subsampler 220 and a Savitsky-Golay smoothing 
routine 222 before taking derivatives and characterizing the data to create a model 
caUbration profile for each ROI being scanned. The model caUbration profiles are 

25 then optionally stored 230 for later optional archiving and restoration 232 under the 
control of the Instrument Control and Data Acquisition Software 202. 

Next, the Instrument Control and Data Acquisition Software 202 is utilized to 
get and set internal parameters 250 in the EPF Fitting module 260. Raw data 262 
from a run-time scan for one ROI performed utilizing the chip is sent firom the 

30 Instrument Control and Data Acquisition Software 202 to the EPF Fitting module 
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260, which qualifies it, queries 264 the EPF Calibration module 210 for the model 
profile 268 for the chip, and then fits the curve utilizing various matrix routines 270, 
iterating the process when necessary. Results 280 from the fitting process are then 
returned to the Instrument Control and Data Acquisition Software 202, which then 

5 provides them to the user. In the preferred embodiment, the EPF Support module 
290 interacts with both the EPF Calibration module 210 and the EPF Fitting module 
260 to provide various testing and common routines. 

Fig. 3 depicts the graphical output of a calibration scan according to the 
present invention. As seen in Fig. 3, the raw calibration scan (prescan) data for ROI 

10 #n are plotted on graph 3 10 as signal strength in A/D counts 312 versus incidence 
angle ^314, Each circle 320 represents an actual scan data point. The step-size, as 
well as its variability, is exaggerated so that it can be seen that the distances between 
steps are slightly irregular. A continuous curve 322 is drawn through the data points 
320 in order to guide the eye, but plays no role in the actual analysis. 

15 In Fig. 3, graph 330 shows the result of smoothing and resampling of the 

initial profile curve for ROI #n, again plotted as signal A/D counts 332 versus 
incidence angle 6 334. The subsan5)led steps represented by sampling points 340 
are all now equidistant apart, allowing the production of smoothed and resampled 
calibration profile 342. Graph 350 shows the corresponding derivative curve 360, 

20 plotted as A/D counts per degree 370 versus incidence angle 6 372. The data points 
380 are the same as sampling points 340 in graph 330. 

Fig. 4 depicts the graphical output of a run-time scan according to the present 
invention. As seen in Fig. 4, data points 402 obtained during a run-time sc^ of ROI 
#n are plotted on graph 410 as A/D counts 412 versus incidence angle ^414. The 

25 angle steps at which data points 402 are obtained may be, and firequently will be, 
irregular. Graph 420 shows the stored calibration profile 430 and derivative 432 
being used to fit ROI #n on the same scale used for graph 410. The profile 430 is 
preferably also from ROI #n, but may be firom another ROI if desired. 

In Fig. 4, graph 440 shows the calibration profile 430 shifted by fiill steps 

30 until it can be fitted to the run-time scan data points 402, resulting in curve yshft 450. 
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Derivative curve 432 is similarly shifted to produce derivative curve 452 and the 
run-time angle points 460 located. Flat baseline model component 462 is also 
shown. The derived empirical profile is shown in graph 470. The raw data points 
402 are shown as squares and the points 474 from the shifted calibration curve yshft 
5 450, here available at fine angle spacing (0.01 deg), are shown as circles. The 
empirical profile is Umited to a selected fitting window that does not include pruned 
raw scan points 480, 482. The raw data points 402 are fitted with a model consistmg 
of shifted profile yshft 450, derivative term 452 (shown sampled at the data point 
separation), and a baseline offset (not shown.) Solid line y 472 is the resulting 

10 empirical profile, essentially a smoothed and shifted version of actual caUbration 
profile 430, which may be fi-om the same ROI or another. In graph 470, the 
residuals are exaggerated for emphasis. 

Calibration Stage. In the preferred embodiment of the present invention, to 
begin the Calibration stage an SPR scan is performed on an SPR chip at relatively 

15 high angular point density, typically 0.05^ nominal point spacing, and over sufficient 
angular range to encompass the fiiU SPR resonance profiles at all ROIs on the chip. 
The full two-dimensional data set obtained is then transmitted to the calibration 
routine along with a vector providing actual measured encoder angles for each angle 
point. It is not necessary that the angle steps be exactly uniform. 

20 Once there, the individual ROI scans are handled individually. Each is 

extracted fi-om the 2-dimensional array and, in the preferred embodiment, is 
prechecked for quality. Each must actually have a resonance dip of adequate depth, 
must have a minimum peak intensity (outside the dip), and must have only one 
major dip. Angles must be monotonically increasing. Once quaUfied as usable, the 

25 resonance is sub sampled to a higher angular density, with rigidly standardized 
angular point spacing, typically 0.010**. In the preferred embodiment, this 
subsampling is carried out by a series of steps as follows: 

I . Raw scans are optionally presmoothed using Savitsky-Golay or any 
other suitable smoothing technique. 
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2. The quasi-randomly spaced points are interpolated using any suitable 
interpolation method, such as cubic spline or Lagrange, in order to determine 
subsampled points on the desired regular angle grid. If desired, this step can be 
avoided or reduced to simple linear interpolation by use of a very fine calibration 

5 step size. 

3 . The ends of the scan are then optionally extrapolated slightly using 
polynomial fits to the end portions of the scan in order to increase the range slightly 
to permit fiirther filtering without loss of scan range. 

4. The interpolated, extended scan is then optionally smoothed for a second 
10 time, using Savitsky-Golay or another suitable smoothing technique. The Savitsky- 

Golay or other smoothing parameters for each stage are adjustable. 

The regularized, subsampled profiles are then optionally examined to 
determine certain properties, which are checked against control limits to establish the 
acceptability of each ROI profile. In the preferred embodiment, the properties 
15 measured include: 

1 . Nominal Resonance Angle, determined by finding the analytic 
minimum of a polynomial fit to the bottom portion of the resonance. The value 
obtained varies according to the polynomial order used and the point selection, but 
the absolute value is not important. Without this property, only angle shifts can be 

20 reported, not actual angle values. 

2. Nominal width, expressed at Full Width at Half Maximum (FWHM) 

3 . Nominal Fractional Depth of resonance 

4. Maximimi Intensity In Profile 

Other properties whose values may be usefiil may also be measured. These results 
25 are then used to estabUsh an overall quality index for each profile. In a typical 

embodiment, the quality ratings are values such as "Excellent", "Good", "OK", 

"Poor but usable", and/or "Bad." Anything similar is, of course, suitable for use in . 

the present invention. 

The subsampled and regularized profiles for each ROI are inserted into a 
30 two-dimensional array, which may optionally be archived for later use, along with 
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the associated property values, particularly including the nominal resonance angle. 
The original angle set used for the profile measurements is generally no longer 
needed and may now be discarded, unless a future use for it has been identified in a 
particular embodiment of the invention. 
5 Finally, derivative vectors are computed for each of the profiles. In the 

preferred embodiment, this is accompUshed using a Savitsky-Golay first derivative 
filter, but other methods know in the art would be suitable. Specific parameters are 
adjustable, but an order of 3 and a Number Of Points of 9 are typically used. The 
derivatives are also combined into a 2-dimensional array for storage in memory. 

10 These are not typically archived, as they can be easily recomputed when archived 
subsampled profile sets are retrieved, but they may be archived if the need arises. 

Fig. 5 is an operational flowchart of the Calibration stage of the preferred 
embodiment of the present invention. Each of the steps is described in more detail in 
the text that follows. A description of the implementation of the preferred 

1 5 embodiment is presented at the end of the Detailed Description section. 

As seen in Fig. 5, in the preferred embodiment the Calibration stage begins 
with an initiaUze and self-test procedure 505. The prescan SPR scan is performed on 
the chip being profiled in order to acquire the caUbration profile data set 510, which 
is then subjected to a preliminary quahty check 520. If desired, the optional chip 

20 qualification check 525 may be performed at this point. If the quahty is not 
acceptable 530, the process is halted 535 and the user is informed of the problems 
discovered. Otherwise, a caUbration profile is extracted 540 and processed 545 for 
the first ROI. If there are other ROIs to be profiled 550, steps 540 and 545 are 
repeated. Otherwise, Calibration set statistics are computed 555 and optionally 

25 displayed to the user. Next, derivatives of the calibration set model profiles are 
computed 560 and the caUbration set is flagged as valid in memory 565. The 
calibration set may next be optionally stored in memory 570. 

Initialize and Self-Test 505. In the preferred embodiment, the detailed 
operation of the EPF caUbration software is controlled by a set of control parameters. 

30 Default values of these parameters are contained within the program and are set by 
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the calibration initialization routine. Optionally, the control software may retrieve 
the current parameter set and/or set new values. The calibration initiaUzation 
procedure also initiates a complete test of the EPF calibration module from end-to- 
end, using test procedures that employ randomized synthetic data. The test results 

5 are then verified for correctness. In the preferred embodiment, the initialization 
process also includes setting a calibration initiahzation flag. 

Acquire Calibration Profile set 510. The SPR signal is measured for each 
ROI as a function of angle of incidence over a range sufficient to encompass all 
resonances on the chip. Angle steps are maintained nominally equal at a 

10 predetermined value, typically 50 mDeg or approximately 5% of the width (FWHM) 
of the resonance. SPR signals are averaged over the detector pixels defining each 
ROI, which comprise typically 200 to 4000 CCD pixels each. ROI shapes may be of 
various shapes, including rectangular, elliptical, or annular. Signals are usually 
expressed in mean A/D (Analog to Digital) counts averaged over the ROI, but any 

15 consistent units, including, but not limited to, photons, photoelectrons, volts, or 
microwatts may be used. Since the detector integration time is fixed for an 
experiment, it does not matter whether light power (microwatts, photons/sec, etc) or 
integrated flux (photons, counts, microjoules, etc) are used. In the preferred 
embodiment, the fiiU data set is contained in a two-dimensional array, with indices 

20 ofROI index and angle index. A separate vector contains the angle values, 

preferably angles measured using an angular encoder. Thus actual measured angles, 
rather than nominal commanded angles, are preferably employed. 

In normal operation, these data are acquired by the instrument control 
software and passed immediately to the main EPF calibration routine. However, the 

25 EPF calibration routine may also be used in a post-run mode for reanalysis of data 
sets, in which case the raw calibration scans will have been stored and retrieved. 

Preliminary Quality Check 520, Routine checks of the appropriateness of the 
data set are made. This step is not essential, but serves to uncover possible problems 
with the data set. Various procedures verify that each ROI label is unique, that the 

30 angle values are monotonically increasing, and that the individual angle step sizes 
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are within spec. One of these routines also estabUshes the precise angle range and 
evenly spaced angle values to be used for the final calibration profile set. 

Chip Qualification Check 525. In one embodiment of the invention, the SPR 
chip itself may optionally be tested against expected quality standards to ensure that 
5 the chip or its preparation is not defective. This optional aspect of the invention is 
described in additional detail below, in conjunction with the description of Fig. 12. 

Extract Calibration Profile Scan 540, In the preferred embodiment, a vector 
extraction routine copies one ROI scan firom the array for detailed processing. 
Process one Calibration Profile (core procedure) 545, This procedure is 
10 described in more detail later, in conjunction with the description of Fig. 6. 

Compute Calibration Set Statistics 555. In the preferred embodiment of the 
invention, certain global statistics are computed, including the mean, minimum and 
maximum values of the following properties: estimated angles of SPR resonance 
minima (Nominal Resonance Angle), approximate FWHM of resonances (FWHM), 
1 5 peak intensity in calibration profile (Maximum Intensity), and fractional depth from 
upper tangent line (Fractional Depth). The global statistics can be retrieved, but are 
not essential in routine operation. 

Compute Derivatives of Calibration Set Model Profiles 560. In the preferred 
embodiment, derivatives of the finished profiles are computed by using the Savitsky- 
20 Golay method, with parameters derivative order and derivative number of points 
typically set to 3 and 9 respectively. 

Optionally store Calibration Set to disc 570. The finished calibration set 
may optionally be archived by the Instrument Control Software for later retrieval and 
use. In the preferred embodiment, derivatives are not archived as part of the finished 
25 caUbration set, but are recalculated when needed. Alternatively, the raw calibration 
scans can be stored and used to recompute the calibration from scratch at a later 
time. 

Fig, 6 is an operational flowchart of the processing of a single ROI 
caUbration profile, shown as step 545 in Fig. 5. Each of the steps is described in 
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more detail in the text that follows. A description of the implementation of the 
preferred embodiment is presented at the end of the Detailed Description section. 

As seen in Fig. 6, the selected ROI calibration scan is checked for suitability 
and key quality features 605. If the scan is found to be unsuitable 610, the ROI is 

5 returned as "Bad" 615. Otherwise, a preliminary smooth is performed 620 on the 
raw profile data, and the resulting curve is subsampled 625 to prescribed angle 
intervals. The ends of the scan are then extrapolated 630 as required to fill out the 
prescribed angle range and a final smooth 635 is performed to generate a member of 
flie finished calibration data set. Next, various properties of the finished calibration 

10 profile are determined 640 and checked against specification limits to determine the 
quality 645 of the finished calibration profile. Finally, the finished profile is stored 
650 to a 2-dimensional array. 

Check single ROI Calibration profile scan for suitability 605. In the 
preferred embodiment, a simple set of rough suitability tests is performed on the 

15 current raw scan. These include, but are not limited to, checking for a sufficient 
number of data poiats, checking the peak and minimum signals in the scan against 
specification limits, testing the ratio of minimum signal to peak signal against the 
maximum limit in order to preliminarily verify an adequate resonance depth, finding 
the angle step at which the minimum signal occurs (presumed to be flie approximate 

20 resonance position), and checking for adequate scan range on either side of this 
minimum. 

In the preferred embodiment, an additional test is made to ensure against 
multiple deep dips that could cause ambiguity in the resonance determination. 
Failure of any of these tests causes the ROI profile property to be flagged "Bad", a 
25 reason why the profile is bad to be determined, and further processing of this ROI to 
be terminated. Warning levels of these quantities will cause the ROI to be flagged 
"Poor", but processing will continue. Failure of any one ROI scan does not prevent 
successful processing of other ROI scans. 

Preliminary Smooth 620. In the preferred embodiment, actual preliminary 
30 smoothing is carried out by the conventional Savitsky-Golay smoothing method. 
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The implementation used in the preferred embodiment, imlike some Savitsky-Golay 
implementations, does not truncate the length of the scan. This is preferably a 
relatively gentle smooth, with parameters 2, 3. This smoothing operation may be 
omitted entirely on clean data without compromising the overall process. Other 
5 suitable smoothing methods known in the art may be similarly utilized. 

Subsample to prescribed angle intervals 625 , In the preferred embodiment, 
the scan is then resampled to the earlier defined angle steps using an interpolation 
procedure. Although Lagrange interpolation may be advantageously used, Cubic 
Spline interpolation performs adequately, and this is employed in the currently 

10 preferred embodiment. Any other method known in the art would also be suitable. . 
The subsampled step intervals are a prescribed subdivision, typically 5, of the 
nominal angle step size (typically SOmDeg). 

Extrapolate ends of scan 630. In some cases, the end points of the final angle 
range for the calibration set are set to be sUghtly beyond the minimum and/or 

15 maximum angles actually measured. In such cases, the end points must be 
extrapolated. The cubic splines determined in the previous step are unsuited for 
extrapolation, so, in the preferred embodiment, a separate polynomial fitting 
procedure is used on one or both ends of the scan, as needed. A polynomial of order 
n is fitted to the final n+1 resampled points determined previously, and this 

20 polynomial is evaluated at the missing angle values (usually no more than one) to 
carry out the extrapolation. Use of this procedure assures that the range of the 
calibration profile is not unnecessarily reduced by truncation due to the inevitable 
mismatch of experimental angle values and the regularized angle steps. Any other 
suitable extrapolation techniques known in the art may also be advantageously 

25 utilized for this step. 

Final Smooth 635 . In the preferred embodiment, the smoothing procedure is 
used again to perform a non-truncating final smooth of the resampled profile, with 
parameters order and number of points typically being 3 and 9. 

Determine Properties of finished cahbration profile 640 . In the preferred 

30 embodiment, tiie detailed properties of the final resampled profile are determined by 
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a separate property determination procedure. Properties determined may include, 
but are not limited to, FWHM, Nominal Resonance Angle, Fractional Depth, and 
Maximum Intensity. Aside from the Nominal Resonance Angle, these properties are 
currently used only for assessing system performance, and for estimating the quality 

5 of the calibration profiles. 

Check Properties against specification limits 645 . In the preferred 
embodiment, the properties Maximum Intensity, Nominal FWHM, and Nominal 
Fractional Depth are checked against specification limits, and the profile quality is 
set accordingly to appropriate quality labels such as "Excellent", "Good", "Poor", or 

10 "Bad". 

Store Finished Calibration Profile 650 . In the preferred embodiment, the 
completed resampled, smoothed and qualified profile is saved to a persistent internal 
data structure. 

The Profile Properties returned by the calibration routines are defmed in 

15 special ways. The defmitions may best be understood with reference to Fig. 7, a 
diagram depicting the definitions of the profile properties returned by epfCal in the 
preferred embodiment of the invention. Calculations are carried out on the 
subsampled profile 702, not on the raw data indicated by the diamond symbols 705, 
First, the minimimi 710 of the curve 702 is identified and is denoted ^dip- In 

20 the preferred embodiment of the mvention, this Nominal Resonance Angle 710 (SPR 
"minimum") is located by a polynomial fit to the bottom region of the subsampled 
profile. Next, the line tangent to the profile 702 from above, called the convex Ud 
715, is determined in order to assist in the estimation of the width and depth of the 
resonance in the presence of a curved baseline. Con5)utations of width and 

25 fractional depth are carried out on the residuals obtained by subtracting the profile 
702 from the convex lid 715. The point 720 shown on the convex lid is directly 
above the SPR minimum 710 and is used to estimate the fractional depth. As shown 
in Fig. 7, the Peak Intensity 725 in the profile 710 is the highest subsampled point 
A line 730, parallel to the convex lid 715, is drawn through the minimum 

30 point 710, and another parallel 735 is drawn midway between the convex lid and the 
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parallel through the minimum. The full width at half maximum 740 (FWHM) of the 
profile is the angular spread between the intersections of the midway parallel 735 
and the side of the resonance profile 702. No attempt is made to interpolate between 
subsampled points, so this value is rounded to the subsampling angle interval. 

5 The firactional depth of the SPR is taken as yi/yo, where yi 745 and yo 750 are 

as shown in the figure. Again, no attempt is made to interpolate between 
subsampled points, so this value is slightly approximated. 
Properties determined are therefore: 
a) Peak Signal in scan 725. 

10 b) Nominal Resonance Angle in degrees 710. In the preferred embodiment, 

this is determined by fitting a polynomial of specified order to the bottom of the 
resonance dip over a specified angle range. Although not fundamentally 
meaningful, this value can be used later to assign self-consistent absolute dip angles 
to fitted scans. 

15 c) The approximate Full Width at Half Maximum (FWHM) of the resonance 

in degrees 740. In the preferred embodiment, to determine the approximate FWHM, 
a procedure effectively passes a line through the resonance minimum that is parallel 
to the convex lid. It then finds the points m the profile on either side of the 
resonance that come closest to midway between the upper and lower lines. The 

20 angular distance between these two profile points is taken as the FWHM. 
Interpolation is not performed. 

d) The Fractional Depth of the resonance. In the preferred embodiment, the 
point on this line directly above the minimum in the resonance is calculated as an 
estimate of the signal intensity Iq in the absence of the resonance. The Fractional 

25 Depth is then determined as the ratio (Iq - Imin ) / Io > where Imin is the intensity at 
the resonance minimum. 

Fig. 8 is a screenshot depicting an example user interface for the 
configuration and initiation of the calibration stage. The same screen is also used for 
the configuration and initiation of the Fitting stage. It must be appreciated tiiat the 

30 screen embodiment depicted in Fig. 8 is provided by way of example only and that 
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additional parameters may be added, parameters may be deleted, or alternate screen 
configurations may be employed without departing from the spirit of the invention. 
The example configuration screen shown in Fig. 8 is used to set up various 
angle scan parameters. As seen in Fig. 8, the following parameters may be specified 
5 and/or verified by the user using the example screen: 

a. Prescan angle range (this is the Calibration profile scan) - "Prescan Start 
Angle" 805, "Prescan End Angle" 810 

b. Prescan angle step size - "Prescan Angle Delta" 815 

c. Scan Extension Factor 820: This is the way the run-time scan range is 
10 specified in the currently preferred embodiment. The average width FWHM of the 

profiles determined in the calibration run is multiplied by this factor to determine the 
Half-range of the run-time scans. For example, if the mean FWHM is 0.8 deg, and 
the Scan Extension Factor is 1 .5, then scans will extend 1 .5*0.8=1 .2 degrees above 
the highest observed ROI resonance, and 1 .2 degrees below the lowest observed ROI 
15 resonance. This scan range may be fixed for the duration of the run, or may be 
adaptive, the scan limits automatically increasing or decreasing as the resonances are 
observed to move about under the influence of refiractive index changes and 
chemical binding. 

d. Step Size Fraction 830: This is a factor that currently determines the run- 
20 time scan angle step size. In tfie preferred embodiment, rather ttian specifying it 

directly, it is calculated as a fi:action of the mean FWHM. It may alternatively be 
specified in degrees, or in any other suitable manner known in the art. 

e. The Bi-Directional checkbox 840, if checked, causes data to be acquired as 
the instrument's Ught source scans both up and down in angle, thereby saving the 

25 time otherwise required to mechanically slew back to the beginning for each scan. 
Use of this mode requires that the Time --of resonance feature of the EPF be 
employed, since otherwise large alternation in the time intervals between resonance 
measurements on successive scans results when the resonances are not exactly at 
midscan. 

30 f. Continuous vs. Stepwise 850: The preferred embodiment of the system 
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supports either or both stepwise start-stop angle scanning or a non-stop measure-on- 
the-fly mode. 

g. Images per Step 860: In stepwise scanning, a variable number of CCD 
exposures (frames) may be measured at each step. Usually one image is obtained 

5 per step, but additional images may be co-added to improve signal-to-noise 

h. Camera Exposure Time 870: The exposure time for each CCD frame is 
chosen to use as much of the CCD pixel well capacity as possible, while avoiding 
any risk of saturation on individual pixels. It depends entirely on the LED Ught 
source intensity, but is typically 25-100 ms. 

10 It should be understood that this and otiier types of conjBguration screens may 

be expanded to include any other suitable calibration and/or fit parameter settings. 

Fig. 9 is a screenshot depicting an example output from the calibration stage. 
It must be appreciated that the screen embodiment depicted in Fig. 9 is provided by 
way of example only and that additional parameters or graphics may be added, 

15 parameters or graphics may be deleted, or alternate screen configurations may be 
employed without departing from the spirit of the invention. In the preferred 
embodiment, this screen is visible while performing the Calibration Profile Scan, 
being updated as each ROI is analyzed. 

As shown in Fig. 9, the button 905 that initiates a calibration scan is shown 

20 on the left. Visible at the top of the screen is a listing of the run method 910, 

selected by clicking on the "Flow Program" tab 911, showing the sequence of fluidic 
control steps which occur during the run and which influence the reaction sequence 
giving rise to the SPR signal. The run temperature (temperature of the chamber 
containing the sample and SPR chip) is shown in box 912. Alternate information 

25 display choices may be had by selection of the "ROI Information" tab 913 or the 
"Result file" tab 914. On the left is an image 915 of the SPR chip itself, brought 
onto the screen by clicking the Image tab 918, showiag the locations of active ROI 
zones 920. Several of the ROIs are highlighted 925; these highlighted ROIs are the 
ones for which angle scans 928 are displayed in the large graph 930 at the lower 

30 right. An alternative to the Image tab 9 1 8 is the ROIs tab 935, in which the locations 
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of the ROIs are displayed graphically, but without an actual monochrome chip 
image. In one embodiment, a similar screen is used to display the results of the run- 
time scan, initiated by clicking the Start Run button 950. 

In the preferred embodiment of the invention, during the calibration scan, 

5 which takes a few seconds, the graph 930 is updated in real time. The actual 
measured intensities, averaged over ROI pixels, are the dot symbols 960 on the plot. 
The lines 928 between them are interpolated for cosmetic purposes using standard 
graphic smoothing functions, such as, but not limited to, splines or bezier curves. 
Note that the ROIs displayed in Fig. 9 exhibit some variation in shape, especially in 

10 the wings, as well as some dispersion in angular position. This is normal. These, and 
the rest of the similar curves for the other ROIs, are the grist for the calibration mill 
In the preferred embodiment, the intensity for each pixel is converted by a 12-bit 
AID converter to an A/D count between a range of 0 and 4095 (2*^). The angle is 
obtained from precision angular encoder readings. 

15 Fig. 10 is a screenshot depicting a different example output from the 

calibration stage, wherein there is a "bad ROI" 41 showing no resonance. It must be 
q>preciated that the screen embodiment depicted in Fig. 10 is provided by way of 
example only and that additional parameters or graphics may be added, parameters 
or graphics may be deleted, or alternate screen configurations may be employed 

20 without departing from the spirit of the invention. 

As shown in Fig 10, the EPF calibration summary 1010 is shown at upper 
right. It is provided when the "ROI Information" tab 1015 is selected. The ROIs tab 
1020 is selected, so the chip is displayed in the graph 1025 at lower left as a series of 
spots 1030 representing the active ROIs. ROIs marked with squares 1035 are 

25 selected for plotting in the chart pane at right. The selected ROI 1040 (ROI 41) is 
"bad", showing no resonance, as can be seen in the EPF Profile Properties window 
1045 and in the resonance curve 1050 in the large graph 1060 at lower right. 

Fig. 1 1 is a screenshot depicting a different example output from the 
calibration stage, showing the details of "good ROI" 190. Again, it must be 

30 appreciated that the screen embodiment depicted in Fig. 1 1 is provided by way of 
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example only and that additional parameters or graphics may be added, parameters 
or graphics may be deleted, or alternate screen configurations may be employed 
without departing from the spirit of the invention. 

As shown in Fig. 1 1, the EPF calibration summary 1 1 10 is shown at upper 

5 right and the EPF Profile Properties 1120 of selected ROI 190 1030 are shown at 
upper left. The ROI tab 1 140 is selected, so the chip is again displayed in the graph 
1 150 at lower left as a series of spots 1 160 representing the active ROIs. Several of 
the ROIs are highlighted; these highhghted ROIs are the ones for which profiles 
1 170 are optionally displayed in the large graph 1 180 at the lower right. 

10 Chip Qualification Chech The mechanisms of the EPF permit a rigorous 

assessment of chip quality to be made from the calibration data, before proceeding 
with the calibration data set generation. SPR chips of particular types or batches 
should have very similar SPR profiles. Any deviations from these profile shapes and 
angular positions are a sign that the chip was incorrectly manufactured or processed, 

1 5 or that it has been damaged. Valuable user time and sample material can be saved 
by identifying defective chips before commencement of the fiill run. 

To implement this feature in the preferred embodiment, representative 
caUbration profiles, averaged over several full chips of known good quality, are 
generated for each type or distinct batch of devices. Although these could be 

20 generated on an ROI-by-ROI basis, in general a single profile representative of the 
whole chip will be used. Expected and allowable variations from this "Golden 
Profile" will also be established. This Golden Profile data, along with the allowable 
variance specifications, can be shipped with the SPR chips from the factory or be 
otherwise distributed to the user, and then stored on the instrument control computer. 

25 Fig. 12 is an operation flowchart for the optional chip quaUfication 

procedure, as shovm in step 525 of Fig. 5. Each of the steps is described in more 
detail in the text that follows. A description of the implementation of the preferred 
embodiment is presented at the end of the Detailed Description section. 

As shown in Fig. 12, the chip qualification procedure is initiated by loading 

30 the correct "golden calibration" 1205 from among the stored golden calibration 
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profiles 1210. Next the fit module is initialized 1220 with the stored chip 
qualification parameter set 1225. A single ROI scan is extracted from the calibration 
scan data 1230, it is processed 1240 as described below in conjunction with Fig. 14, 
and the results are analyzed for conformance with the golden profile 1250. If the 
5 results for that ROI fail the test 1260, the count of bad ROIs is incremented 1265 and 
the number of bad ROIs is optionally stored for reporting to the user. If additional 
ROIs remain to be tested 1270, another single ROI scan is extracted from the 
calibration scan data 1230, and the process repeats until all ROIs have been 
processed. At flie end of testing, if the number of bad ROIs found exceed the 
10 predetermined threshold 1280, the discrepancies are displayed to the user 1285 
and/or the chip is rejected 1290. Otherwise, the chip qualification procedure returns 
an indication that the tested chip is within acceptable limits 1295. 

The following steps are taken to implement the Chip Quahfication Check 

(CQC): 

15 Load Golden Calibration . The stored Golden Calibration archive 

corresponding to the chip type or batch number is loaded and made to be the current 
calibration set. It will typically consist of a single profile, rather than a large set of 
profiles for each ROI. 

Initialize epfFit with Special Parameter Set. The epfFit parameters used for 

20 the CQC are generally different from those used in normal analysis operations. In 
particular, the specification Umits for the various parameters will vary, some being 
tighter than normal, and others being looser. In particular, the allowed rms 
Residuals will be larger, since the individual ROI scans are being compared to a 
single standard profile rather than to individuahzed ROI calibrations measured on 

25 the same ROI of the same chip. 

Extract single ROI scan from Prescan Data . The same calibration data set 
that is to be used for the calibration stage is used for the CQC. One ROI at a time is 
extracted from that set and submitted to the specially configured epfFit procedure. 
Note that ttie data point spacing is that of the caUbration scan, and is therefore 

30 smaller than the spacing typically used for run-time scans. 
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Process one ROI Scan , This is identical to the process used for fitting 
ordinary run-time data, with two exceptions: 

a) The epfFit parameters are modified, as discussed above, 

b) The Calibration Profile to be associated with each ROI is different. 
5 Instead of using the recently measured Calibration Profile for the ROI in question, 

the Golden Profile is used for all ROIs. In the preferred embodiment, the ROI to 
Use is one of the calling parameters in epfFit. In this case, the value is "Golden 
Profile", or whatever label has been used for the Golden Profile. 

Analyze ROI Results for Conformance with Golden Profile . The normal 

10 fitting process, with the modified specification limits, generates a number of quality 
checks in itself The rms Residual value is the most important of these, and is 
therefore the primary test criterion, as it detects any significant change of shape of 
the SPR profiles from the expected standard shape, including changes in FWHM. 
As manufacturing and processing consistency changes, the allowed variances on this 

1 5 and other parameters may be modified. Appropriate limits on the key parameters are 
typically established in the light of field experience, in order to avoid excessive false 
rejections. 

However, other parameters, not usually of direct concern in normal fitting, 
should also be checked. These include: 

20 a) The SPR angle shift. Since the measurement is being made prior to 

the run, SPR angles should vary fi*om the golden standard because of two effects: (i) 
the refractive mdex of the buffer solution in use, and (ii) the surface treatments or 
ligands applied to the chip by the user. Item (i) can be explicitly allowed for, given a 
definition of the buffer in use. Item (ii) cannot in general be quantitatively taken into 

25 account, except that the effect of (ii) is always to increase the SPR angle. Therefore, 
after taking (i) into account, any negative angular shift (beyond a specified tolerance) 
is a sign of trouble. 

b) Baseline Fit CoeflBcient. If the profile shape, including resonance 
depth, is close to that of the golden profile, then the Baseline component in the fit 
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should be very small. Therefore, placing limits on this fit result is equivalent to 
ensuring that resonance depth is within allowed tolerances. 

c) Profile Fit Coefficient. This value should be of order unity, but may 
be allowed to dip below unity in order to account for reduced signal intensity on the 
5 target instniment. The golden profiles will have been normalized to full intensity, so 
this fit coefficient should be allowed to exceed unity by only a small amoimt. 

When an ROI fails one or more of these tests, it is flagged as "Bad", and the 
cause of the discrepancy logged. A count is kept of the number of such out-of-spec 
ROIs. 

10 Test for Too Many Bad ROIs. When all the ROIs have been analyzed, the 

number of Bad ROIs is tested against a predetermined limit, which may be zero. If 
the number exceeds this limit, the chip is rejected - or, at a minimum, the user is 
informed and given the option of terminating the run before committing sample 
solution. If desired, the user may also be informed of the specific ROIs that are 

15 "Bad", in order to allow the user to decide if the affected ROIs are of particular 
interest in the planned experiment. 

Fitting Stage. In the preferred embodiment, when resonances obtained in 
the course of an experiment are to be measured, the epfFit routine is called. In this 
case, just one profile, corresponding to a single ROI, is submitted at a time. Three 

20 vectors are submitted: the SPR intensities, the encoder angles, and the clock times 
for each data point, along with the ROI ID to be used in fitting. 

The angles at which the profile is sampled are generally not the same as those 
used during caUbration, nor are they necessarily spaced similarly. Typically, much 
sparser angle spacing is used in order to speed up measurements. Wider angular 

25 spacing has little effect on angle determination noise until a minimum number of 
points per FWHM of the resonance, typically three or four, is reached, whereupon 
coarser spacing definitely degrades angle accuracy. Shot noise (and other amplitude 
noise) in the data does reduce accuracy, however, and, to the extent that coarser data 
point spacing results in fewer total photons going into the measurement, there is a 

30 deleterious effect. 
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It is possible to obtain valid fits, albeit with reduced precision and accuracy, 
with partial scans (i.e., incomplete resonance profiles). This is useful in extending 
the dynamic range of effective refractive index (RI) values that can be measured. 

The fitting process is more complex than the calibration process outlined 
5 above. The essential steps are: 

1. Data vectors are checked for consistent array indices. Angles are 
checked to make sure that the scan is monotonic in angle, although either increasing 
or decreasing angle is acceptable. 

2. The ROI designation is checked to be sure that a calibration exists for it. 
10 It is not essential that the calibration profile used for fitting belong to the same ROI, 

although normally it does. If for some reason this ROI was never caUbrated, or the 
calibration failed, then another similar ROI calibration profile may be used instead 
with reasonable success. This is because the fit results are not extremely sensitive to 
profile shape, and because it is not the absolute SPR angle, but rather shifts of angle, 
1 5 that are of primary interest during the course of an experiment. 

3. If the angles are in decreasing order, the data vectors are inverted for 
convenience in comparing wifli the calibration profiles, which are in strictly 
increasing angle order. 

4. The rough location of the SPR minimum is determined, using a judicious 
20 combination of locating the lowest point in the scan and fitting a small number of 

nearby points with a low order polynomial (quadratic). If no clear minimum can be 
found, the procedure fails, or, in some cases, the lowest point (if at one end of the 
scan) is used as a first guess. 

5. The correct subsampled calibration profile and derivative, previously 
25 computed in the calibration stage, are retrieved. 

6. The calibration profile and derivative are shifted in full subsample steps 
(typically 0.0 r) to bring the calibration profile nominal minimum into coincidence 
with the rough SPR location determined in step (4) above. This is merely an index 
pointer calculation. 
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7. The overlap portion of the SPR scan to be fitted and the calibration 
profile is identified. In general, the overlap is shorter than either alone, but 
sometimes one is contained within the other. Only this overlap portion is used. 

8. Signal and derivative values at the actual SPR scan angle values are 
5 derived from the shifted calibration profile. In the preferred embodiment, this is 

accomplished with Lagrange interpolation, but any other suitable method known in 
the art may be xised. This is accurate because the angle point spacing in the 
calibration profiles is very small. 

9. The pruned portion of the SPR scan is fitted. In the preferred 

10 embodiment, this is accomplished using classical linear least squares fitting, using as 
model components the following: 

a. The calibration profile 

b. The calibration profile derivative, 

c. A constant additive offset. 

15 10. The fractional angular shift of the apparent SPR minimum from the 

shifted calibration profile is determined as the quotient of two of the fit coefficients 
obtained in step (9) above, namely the ratio of the calibration derivative fit 
coefficient to the calibration profile coefficient. The latter is normally 
approximately unity, whereas the former is typically very small. This fractional shift 

20 is added to the discrete shift applied in step (6) above, in order to provide an estimate 
of the actual total SPR shift relative to the SPR location at calibration. 

11. If the fractional shift obtained in step (10) exceeds half a subsample step, 
then steps (6)-(10) are repeated (iterated) until the fractional step is equal to or less 
than half a subsample step. Other termination criteria may be used successfully, 

25 although the described criteria are currently preferred. Special means are provided 
to prevent endless jumping to and fro when the fractional shift is very close to half a 
step. 

12. At this stage, the SPR location has been fairly well established. The user 
may optionally specify, however, that a limited region of the SPR scan near the 

30 resonance be used for the final fit in order to reduce errors due to artifacts in the 
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wings of the profile. When this "Sweet Zone'' fitting is used, the SPR scan is further 
pruned to the specified reduced width, centered as much as possible on the resonance 
location, and the process from step (6) to step (1 1) is iterated again. The reason for 
not pruning to the Sweet Zone the first time around is to reduce the Ukelihood of 

5 grossly incorrect fits by using as much data as possible. Once the correct alignment 
of calibration and unknown profiles has been fairly accurately determined, further 
refinement using pruned scans is safe, even in pathological cases such as highly 
truncated resonances. However, higher execution speeds may be obtained by 
skipping directly to the pruned fits. 

10 13. Various quality checks are performed. In particular, at least three points 

must be remaining after overlap and any further pruning in order to perform valid 
fits. The fitting residuals are checked against both lower and upper limits. 

14. In addition to the angle shift determined above, an estimate of absolute 
resonance angle is made by adding this shift to the previously determined 

15 Calibration Nominal SPR Minimum angle. 

15. The clock time at which the minimum occurred in the scan is estimated 
by interpolation between the times at which the scan data points on either side of the 
found minimum were measured. This Time of Minimum is highly useful for more 
accurately assigning the SPR angle obtained to a specific clock time for use in 

20 kinetic analysis. 

16. Results of the fit are returned to the calling program, along with various 
details such as fit coefficients and a quality estimate based on comparing these 
parameters with control limits. 

Fig. 13 is an operational flowchart of the fitting stage of the present 
25 invention. Each of the steps is described in more detail in the text that follows. A 
description of the implementation of the preferred embodiment is presented at the 
end of the Detailed Description section. 

As shown in Fig. 13, in the preferred embodiment the fit stage is begun with 
an initialize and self-test procedure 1310. Then the SPR run is performed and the 
30 run-time angle scan at time t is acquired 1320. Next, a single ROT scan is extracted 
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1330 and processed 1340 for the first ROI. If there are other ROIs to be profiled 
1350, steps 1330 and 1340 are repeated. Otherwise, the graphical display is updated 
1360 and the fit results are stored in memory 1370. If another angle scan is desired 
1380, the procedure returns to step 1320; otherwise, the fit stage completes. 

5 Initialize and Self-Test 1310. In the preferred embodiment, the detailed 

operation of the EPF fit software is controlled by a set of control parameters. 
Default values of these parameters are contained within the program and set by the 
fit initialization routine. Optionally, the control software may retrieve the current 
parameter set and/or set new values. The initialization routine also initiates a 

10 complete test of the EPF Fit module firom end-to-end using various test procedures, 
which in turn employ randomized synthetic data generated by a procedure from the 
calibration module. The test results are verified for correctness. In the preferred 
embodiment, the initiaUzation process also includes setting the fit initialization flag. 
Perform Run, After acquiring a calibration set, a run is generally performed 

15 in which a sequence of angle scans is made in order to track the time evolution of the 
SPR angle responses of the several (typically 100 - 400 or more) ROIs on the chip as 
various buffer and sample solutions are flowed over it. The resulting SPR responses 
will later be analyzed to determine the events and parameters of interest, such as 
binding affinities, kinetic rate constants, etc. 

20 Acquire Runtime Angle Scan at time 1 1320. For each time point during an 

SPR run, which may last firom a minute or two to several hours or more (but is 
typically 20-200 minutes, depending on the time scales of the chemical reaction rates 
involved), an angle scan is performed, followed by EPF fits to determine the 
resonance angles for each of the ROIs. Fitting is typically overlapped in time with 

25 data acquisition of the next time point. 

The SPR signal is measured for each ROI as a fimction of angle of incidence 
over a range sufficient to include all resonances on the chip. This range may be 
similar to, greater than, or less than the angle range used for the calibration step. In 
the preferred embodiment, it may be changed during the run as needed in order to 

30 ensure inclusion of all the ROI resonances as resonances shift in the course of the 
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run. The angle steps are maintained nominally equal at a predetermined value, 
typically 50 to 200 mDeg or approximately 5-20% of the width (FWHM) of the 
resonance. Typically, these angle steps are much larger than the steps needed for the 
calibration procedure. 
5 As in the calibration stage, SPR signals are averaged over the detector pixels 

defining each ROI, which comprise typically 200 to 4000 pixels each. ROIs may be 
of any of the various shapes known in the art, including rectangular, elUptical, or 
annular. As discussed, signals are usually expressed in A/D counts, but any 
consistent units may be used. In the preferred embodunent, the full signal intensity 
10 data set is contained in a two-dimensional array with indices of ROI index and angle 
index. 

A separate vector contains the angle values, preferably angles measured 
using an angular encoder. Thus actual measured angles, rather than nominal 
commanded angles, are preferably employed. Another separate vector, not needed 

15 for the calibration data set, contains the mean actual clock times at which the 
iadividual angle data points (each consisting of one or more CCD detector frames) 
were obtained. , a single CCD exposure is used at each angle value, but multiple 
frames may be averaged or co-added for enhanced signal/noise. 

In the preferred embodiment, during normal operation these data are acquired 

20 by the Instrument Control Software and are passed to the main EPF fit routine 
immediately following each angle scan. However, the EPF may also be used in a 
post-run mode for reanalysis of data sets, in which case the measured data will have 
been stored to disc and retrieved. 

Extract single ROI scan 1330. Since the epfPit routine of the preferred 

25 embodiment processes only one ROI at a time (unlike epfCal, which handles the 
whole set of ROIs in a single call), the data for each ROI must be separated out by 
the Instrument Control Software. For each ROI, the ROI index and/or label, a vector 
of pixel-averaged signal intensities, a vector of encoder angles, and a vector of frame 
measurement times is prepared and employed for the fit routine call. Note that the 
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last two vectors remain the same across all ROIs for each scan, but will change from 
scan to scan. 

Process one Run-Time ROIScan 1340. This core procedure is described in 
more detail later, in conjunction with the description of Fig. 14. 

5 Update Graphical Display 1360, In the preferred embodiment, the 

Instrument Control Software generally displays SPR response (angle) vs. tune 
curves for several or all ROIs in real time so that users can observe the chemical 
binding curves developing, identify problems, and possibly determine whether to 
terminate the run early or to extend it. 

10 Store Fit Results to disc 1370, Fit results, including quality measures and 

error codes, are optionally stored to disc for later detailed analysis. In the preferred 
embodiment, for each ROI run, results consist of a vector containing, at each time 
data point, a Results structure containmg the following: epfFit Results, ROI label, 
average step spacing in scan (in degrees), scan inversion flag (True if scan was 

15 downward in angle and had to be inverted), angular shift A6 from calibration profile 
to this scan, ^spr based on nominal calibration profile SPR angle, run time at which 
resonance minimum was observed, number of angle scan points actually used in fit, 
RMS residual (SyRMs), amplitude of fit (should be of order unity), Baseline of fit, 
extrapolation flag (True if result was extrapolated due to incomplete profile), scan 

20 pruned flag (indicates that less than fiill scan lengfli fitted (usual case)), a Quality of 
Result parameter (typical possible values being "Excellent", "Good", "Poor", and 
"Bad"), a description of the primary error condition, if any, and the number of non- 
fatal warning conditions encountered. 

Fig. 14 is an operational flowchart of the processing of one ROI scan, as 

25 shown in step 1340 of Fig. 13. Each of the steps is described in more detail in the 
text that follows. A description of the implementation of the preferred embodiment 
is presented at the end of the Detailed Description section. 

As shown m Fig. 14, initially the availability of a cahbration profile for the 
ROI and the self-consistency of the extracted data vectors are checked 1405, the scan 

30 indexing is checked 1410, and the SPR dip angle is roughly estimated 1415 to make 
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sure it is within the scan range. If any of these tests fail 1420, the processing is 
stopped and an error is returned 1425. Otherwise, the calibration model profile for 
the ROI is fetched 1430 and the SPR angle is determined by fitting to the calibration 
model 1435. If a Sweet Zone fit is requested 1440, the scan is truncated 1445 and 
5 the SPR angle is recalculated by fitting flie truncated scan to the calibration model 
1450. Finally, the time of occurrence of tfie SPR minimum is estimated 1455. 

Check for Calibration Availability 1405. In the preferred embodiment, a 
consistency check procedure is used to verify that the signal intensity, scan angle, 
and time of measurement vectors are all consistently indexed. Next, it is verified 

10 that a current calibration set exists and corresponds to the Chip ID of the current 
data, and that a calibration profile for the specified ROI exists and is valid. If there 
is a problem, appropriate error flags for this ROI are set. 

Check Scan Indexing 1410 . In the preferred embodiment, if the sense of the 
scan is reversed (i.e., from large to small angles of incidence), then a working copy 

15 of the three data vectors with the correct sense is generated. This may happen on 
alternate scans when using bi-directional scanning, which speeds up data acquisition. 
This routine also checks that the angle steps are monotonic and sufficiently 
consistent in magnitude, that there are a sufficient number of angular data points, 
that a sufficient angle range is covered by the scan, and that the time of measurement 

20 values are monotonic in time. 

Roughly Estimate SPR Angle 1415. Initially, a first estimate of the angular 
location of the resonance is generated. First, the minimum intensity point in the scan 
is located. Depending on the settings of certain control parameters, the minimum 
may be required to fall within the scan and at least a certain specified distance from 

25 the ends, or, alternatively, extrapolation may be permitted and the resonance position 
may be allowed to fall outside the scan range. If the minimum is not too close to 
either end of the scan, according to these settings, then the scan is qualified. If the 
minimum is at least one data point from each end of the scan, then a quadratic fit is 
used to estimate an interpolated angular minimum. If not, then the minimum itself is 

30 used as flie rough estimate. This roughly estimated SPR angle is used as a starting 
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point to facilitate the fitting process below. Other methods known in the art for 
roughly estimating the SPR angle would also be suitable. 

Accurate Determination of SPR Angle 1435. 1450 . This core fitting 
procedure is described in more detail later, in conjunction with the description of 

5 Fig. 15. First, the appropriate calibration model profile and derivative vectors are 
fetched from memory. Next, the fit calculation ^Determine SPR Angle by fitting to 
Calibration Model', described in detail conjunction with Fig. 15, is performed once, 
using the full scan range. Then, if a Sweet Zone fit has been requested, the scan 
range is pruned to the specified Sweet Zone range on either side of the first 

10 calculated SPR angle, and the core fitting procedure is repeated to fmd the final 
result. This iteration is performed because the initial fixU range fit is more robust 
than the limited range Sweet Zone fit, and is thus better able to find the approximate 
resonance angle without potentially being "fooled" by anomalies in the data or by a 
highly inaccurate initial estimate. On the other hand, the Sweet Zone fit may be 

15 specified to further refine and optimize the fit, since it is less influenced by possible 
shape artifacts in the far wings of the scan data. 

Truncation to the Sweet Zone 1445 . In the preferred embodiment, truncation 
to the Sweet Zone is carried out in two steps. First, a procedure determines whether 
the scan is akeady too short on one end or the other to encompass the putative Sweet 

20 Zone. If so, the angular limits of the Sweet Zone are shifted away from the short end 
if possible, in order to maintain the full angular width of the Sweet Zone. In the 
second step, another procedure carries out the actual truncation. 

Estimate Time of Minimum 1455 . Knowing the fractional data point at 
which the SPR resonance was found, in the preferred embodiment the time of 

25 measurement is interpolated in order to determine the time m the run at which this 
ROI's SPR minimum occurred. It is not sufficient to simply use the start time or the 
midpoint of the angle scan, since some ROIs may have minima near the beginning of 
the scan and others near the end. Accurate times are important in later analysis of 
kinetic binding curves. 
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Core Fitting Procedure. Fig. 15 is an operational flowchart of the 
determination of the SPR angle utilizing an empirical profile fit, as discussed in steps 
1435 and 1450 of Fig. 14. Each of the steps is described in more detail in the text 
that follows. A description of the implementation of the preferred embodiment is 
5 presented at the end of the Detailed Description section. 

As shown in Fig, 15, the estimated offset of the run-time SPR resonance 
angle from the calibration scan for the same ROI is computed 1 5 10, and the scan is 
repnmed to the calibration profile limits 1520. Next, the cahbration profile is 
interpolated and the derivative model vectors are calculated at the estimated offset 

10 1530. A fit is performed on the repnmed scan using the interpolated calibration 
profile 1540 and the residual angular shift from the estimated shift is calculated from 
the fit coefficients 1550. An improved estimate of the SPR resonance angle is 
calculated 1560, as are the fit residuals 1570. While in the preferred embodiment, 
this fit loop is iterated until the magnitude of the angle shift is less than half the 

15 subsample angle step spacing in the calibration profile 1580, any suitable 

convergence criterion may be advantageously employed. Iteration will take place 
until the value of the resonance shift converges to the predetermined convergence 
criterion. 

In the preferred embodiment, the basic SPR angle determination is carried 
20 out by a procedure, which is called twice. On entry, this procedure is supplied with 
the calibration model profile and derivative, as well as an initial estimate A'spr of the 
SPR resonance angle. On exit, it returns a Results structure. 

A9Rg T- g^g PR-gr Arl510. The estimated offset AGest of the SPR 
resonance from its location in the calibration scan is computed in order to allow 
25 appropriate shiftmg of the model components. 

Reprune Scan to Cal Profile Limits 1520. In general, the current angle scan 
may extend beyond the range covered by the cahbration model on one or both ends. 
Accordingly, the scan is temporarily truncated to a range that does lie within the 
Calibration range, as determined using the current estimated angle offset, AOest- 
30 Note that points temporarily pruned off may be restored on later iterations. 
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Calculate Interpolated Cal Profile 1530. Calculate Interpolated Cal Profile 
1530 is one of the more delicate parts of the fit algorithm. Based upon the current 
estimated angle offset, A9est , the subsampled calibration model profile and its 
derivative are shifted an integer number of subsampled angle steps (which are 
5 typically 10 mDeg) in order to bring them into the closest possible alignment with 
the estimated position of the scan being fitted. Then the model profile and derivative 
are interpolated to the actual irregularly spaced measured angle values of the 
repnmed scan being fitted using Lagrange interpolation. For end points, 2^^ order 
Lagrange interpolation is used, while for interior points in the pruned scan, 4* order 

10 Lagrange interpolation is employed. At this point, resampled calibration profile and 
derivative models have been created at the estimated shift of the scan being fitted. 
While Lagrange interpolation is used in the preferred embodiment, many other 
methods known in the art would be suitable. 

Perform Least Squares Fit of Repruned Scan 1540 . In the preferred 

15 embodiment, the pruned scan is subjected to a classical linear least squares fit using 
three components: 

1 . The resampled calibration profile 

2. The resampled calibration derivative 

3. A constant baseline offset 

20 The third component is not strictly necessary, but is recommended to 

accommodate baseUne system drift and other instrumental effects. Additional 
components can be included to accommodate intensity drift and other instrumental 
artifacts, although some may potentially degrade rather than improve estimates of 
the angle shift. The fit model is: 

25 ypiT = Co ypROFiLE Ci yDERiv + C2 

Calculate 86 from Fit Coefficients 1550. The residual angular shift 5^ from 
the estimated angle used to resample the calibration profile is proportional to the fit 
coefficient for the derivative component. More specifically, 
65 = -Ci/Co 
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where Co and Ci are the first and second model component fit coefficients 
respectively. 

Update Angle Estimate 1560 . The improved estimate of the SPR resonance 
angle is calculated as: 

5/i> NEW ^ >)s I jjA 
SPR *" SPR"^' (>^' 

Calculate Fit Residual 1570 . The residuals are computed as 6yi = Yobserved i 
- Yfit i , and the nns residual value is conq)uted as 
5yRMS =[ E5yi']/DoF 
where the number of degrees of fi-eedom DoF = Npoints - Ncomponents , 
1 0 with Ncomponents ^ 3 . 

|56|< A9rAr/2? 1580. In the preferred embodiment, the iteration of the fit 
loop is terminated when the magnitude of the angle shift bd is less than half the 
subsample angle step spacing in the calibration profile, indicating that the integer 
poiQt shift was as close as possible to the true shift. In most cases, only one pass 
1 5 through the loop is needed, but if the initial guess was poor, one or more further 
iterations may be required. 

The loop termination criteria also include the possibility of oscillation, in 
case the fitted angle value is very nearly halfway between two subsampled 
calibration profile points. It is then possible that the value of \bd\ might stay slightly 
20 above 0.5 AOcal as the integer offset oscillates between two adjacent values on 
successive iterations. In the preferred embodiment, if the algorithm becomes 
trapped ia such a loop, it exits afl:er several iterations, and the lowest observed value 
of is used. Note that although this is the preferred termination criterion, other 
less restrictive termination criteria, such as |5e|< AOcal? can be used with reasonable 
25 success. 

Generate Results 1590. ^spr = spr. The fmal value of the fitted angle 
estimate, on the final pass through the angle fit procedure, is taken as the result. The 
corresponding rms residual is also reported. The Results structure for each ROI is 
created. It will later optionally be stored to disc. 
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Fig. 16 is a screenshot depicting an example output from the fit stage. It must 
be appreciated that the screen embodiment depicted in Fig. 16 is provided by way of 
example only and that additional parameters or graphics may be added, parameters 
or graphics may be deleted, or alternate screen configurations may be employed 

5 without departing fi-om the spirit of the invention. This screen is similar to the 
screens depicted in Figs. 9 - 1 1, but is shown during an actual run, which was 
initiated using the Start Run button 1605 at the left. 

As shown in Fig. 16, the ROIs Tab 1610, rather than Image tab 1615, is 
selected. The various ROIs 1620 shown on the ROIs map 1622 have been spotted 

10 with various ligands. Visible at the top ofthe screen are the EPF profile properties 
for ROI 190 and the EPF prescan summary 1625, selected by clicking on the "ROI 
Information" tab 1630. The run temperature is shown in box 1635. In this case, the 
large graph at the lower right 1640 displays not Intensity vs. Angle, but rather Angle 
Shift vs. Time as the experiment progresses. Again, a subset ofthe ROI curves 1645 

15 is displayed, corresponding to the highlighted ROIs 1650 on the ROIs map 1622. 
Individual data points on this plot each correspond to a complete angle scan, firom 
which EPF calculates a resonance angle position. In this case, new points are plotted 
approximately every 6 seconds, so that the points are too closely spaced to be 
resolved on this scale. Because of this, symbols are not plotted; only the connecting 

20 line segments are displayed. 

During the run displayed, the various ROI plots are autozeroed at the 
beginning of the run, with only the deviations firom the initial angles being 
displayed. At 300 sec, a sample is introduced and binding to the surface commences, 
to different degrees on the various ROIs. At t=900 sec, the sample flow is turned 

25 off, the buffer flow is reinstated, and the bound sample begins to slowly dissociate 
from the surface. 

These sorts of kinetic curves (association phase and dissociation phase) are 
what are fiirther analyzed, in postprocessing, to generate the desired measurements 
of kinetic bmding constants. In the preferred embodiment, this analysis is not done 
30 on the instrument, using the instrument control software, but instead the fit data are 
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exported to other specialized tools. Displays of the finished data can be readily 
generated using these other tools. 

Alternate Configurations. It is feasible to use tabulated theoretical model 
functions in the EPF formahsm, as if they had been measured. In other words, the 

5 calibration scans may be performed as in EPF, theoretical response curves may be 
fitted to these measured curves, and then the fitted curves may be used as the EPF 
model functions. In effect, this provides an alternative means of performing the 
smoothing operation in deriving EPF models. If the experimental profiles are clean 
enough to be accurately modeled by the theory, this works well - but probably not 

10 significantly better than the EPF approach. The computational overhead might 
possibly be acceptable, since the nonlinear fits need be done only once per run, per 
ROI. There is a large increase in sofbsrare complexity, however, with little or no 
benefit. In addition, Explicit Functions, if suitable ones can be identified, may be 
used as stand-ins for empirical profiles, and work about as well as empirical profiles 

15 in the EPF approach. 

Although the above description of the preferred embodiment refers to an 
angle scaimed optical resonance system, as previously discussed, the same 
procedures may be applied to a wavelength-scarmed system. In such case, 
references to "angles" and "angle steps" are replaced by **wavelength" and 

20 "wavelength steps". In addition, references to "encoder angles" are replaced by 
"calibrated wavelength values". It will also be appreciated that the nimierical values 
and units of Wavelength Step sizes are different, such as "0.5 nm" instead of 
"O.OSdeg". The resonance shift is therefore a wavelength shift instead of an angle 
shift. Similarly, the EPF system of the present invention can be applied to data 

25 obtained fi:om predispersion or correlated scanning, with appropriate changes being 
made to the referenced parameters. Further, while use of caUbration profile 
derivatives to determine the resonance shift is currently the preferred embodiment of 
the present invention, other methods known in the art would also be suitable and are 
contemplated as being within the scope of the present invention. 

30 Error Reporting. The preferred embodiment of the invention also includes 
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an error-reporting function. Error reporting may be handled through any of the 
many methods known in the art. Errors may be reported during either the calibration 
or the fitting stage, or both. In one embodiment, an error severity level is specified 
above which an error will be logged. Optionally, a number of errors detected during 
5 the calibration stage that is above a specified threshold may result in rejection of the 
calibration scan or of the chip itself. 

As an example, Fig. 17 is a screenshot depicting an example embodiment of 
a local error log, showing reports of bad fits to particular ROIs. It must be 
appreciated that the screen embodiment dq)icted in Fig. 17 is provided by way of 

10 example only and that additional parameters may be added, parameters may be 
deleted, or alternate screen configurations may be employed without departing fi-om 
the spirit of the invention. As seen in Fig. 17, in this embodiment the error log 
reports the date/time 1710 of error,' the application 1720 generating the error, the 
error severity 1730, the user 1740, the module 1750 generating the error, an indicator 

15 1 760 showing whether the error was caused by bad data received or by some other 
type of event, and a description 1770 of the error. 

Fig. 18 is a screenshot depicting an example embodiment of an error log on a 
remote networked supervisory computer, showing reports of bad fits to particular 
ROIs. It must be appreciated that the screen embodiment depicted in Fig. 1 8 is 

20 provided by way of example only and that additional parameters may be added, 
parameters may be deleted, or alternate screen configurations may be employed 
without departing ft-om the spirit of the invention. As seen in Fig. 18, the error log 
reports the time 1810 of the error, the facility 1820 generating the error, the error 
priority 1830, the apparent error source 1840, the real error source 1850, and the 

25 logged error message 1860. 

Implementation of the Preferred Embodiment. While the present 
invention may be implemented by any suitable method, mechanism, or combination 
of methods and/or mechanisms known in the art, it is preferably implemented in 
software. In the preferred embodiment, the software source code is implemented 

30 using Visual Basic or C++, but any suitable progranrniing language or tool of 
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software implementation known in the art is also within the scope of the invention. 
Similarly, the software maybe run utilizing any suitable operating system, compiler, 
interpreter, application program, or other such device known in the art. In the 
preferred embodiment of the invention, the actual product shipped consists of object 
5 code and associated Dynamic Link Libraries (DLLs). 

Exposed Methods for Calibration. In the preferred embodiment, many 
procedures are Boolean Functions returning True on failure, except as indicated. An 
Error Code of type Long is also generally returned. Many procedures are optional; 
the two that are essential are the calibration initialization procedure and the main 
10 calibration module. The calibration initialization procedure must be called before 
any other action. It initializes tfie calibration module, sets default values for the 
calibration parameters, and performs an end-to-end self test of the calibration 
module. 

The main calibration procedure performs the calibration, storing the results in 

15 its own memory. It does not fail unless a fatal error is encountered. Thus, one or 
more ROIs may fail without causing the whole procedure to collapse. The number 
of ROIs for which a problem was encountered is returned. The problem may range 
from outright failure to a poor rating based on various quality measures. The calling 
program should check whether the number of ROIs for which a problem was 

20 encountered is >0, and, if so, investigate. 

Exposed Methods for Fitting. In the preferred embodiment, only one ROI is 
fitted per call. This allows flexibility, in that other work can be done in between 
ROI fits, such as plotting results incrementally. In addition, this allows more 
detailed error and problem reporting on a per ROI rather than global basis. Fitting is 

25 more complex and tricky than calibration, but there is only one basic method 

involved. Again, most procedures are Boolean Functions returning True on failure. 
Only two procedures are mandatory, the EPF initialization procedure and the main 
fit module. The fit initiahzation procedure must be called before any other action. It 
initializes the fit module and sets default values for the fit parameters. The main fit 

30 routine is a Boolean Function, returning True if any trouble was encountered. This 
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does not necessarily mean that the result is useless; results must instead be inspected 
for details. 

In the preferred embodiment, time is measured in seconds from the initiation 
of a run. In some cases, the times are actually stored as minutes for easy human 

5 interpretation. For example, the affmity plots now have the lower time scale in 
minutes. The measurement time values received by the epfFit routine need not be 
highly accurate. They are used to provide more accurate estimates of the time axis 
for dynamic situations such as kinetic analysis. Probably values accurate to 0,1 sec 
are adequate, although somewhat better values may be useful. These times are 

10 normally reported in seconds, according to some system clock that runs continuously 
during a kinetic run and will be used to plot and analyze kinetic data. Time of day in 
seconds would also be acceptable, except that it will "reset'' at midnight and could 
therefore corrupt a run taken at that time. 

The apparatus and method of the present invention, therefore, provide a new 

1 5 mechanism for quantitation of surface-binding optical resonance curves that has 
relatively low computational complexity, requires a minimal number of scan data 
points, is less sensitive to shot noise, accommodates the use of long scan times, 
provides increased scan speeds, provides results for a sample quickly, accommodates 
scans over a greater number of ROIs, and provides useful results even when only 

20 partial curves are available. Each of the various embodiments described above may 
be combined with other described embodiments in order to provide multiple 
features. Furthermore, while the foregoing describes a number of separate 
embodiments of the apparatus and method of the present invention, what has been 
described herein is merely illustrative of the application of the principles of the 

25 present invention. Other arrangements, methods, , and substitutions by one of 
ordinary skill in the art are therefore also considered to be within the scope of the 
present invention, which is not to be limited except by the claims that follow. 
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CLAIMS 

What is claimed is: 

1 . A method for quantitation of surface-binding optical resonance profiles 
comprising, in combination, the steps of: 

obtaining at least one calibration result from a calibration scan of at least 
one Region of Interest; 

generating, from at least one calibration result, a calibration profile for at 
least one scanned Region of Interest; 

obtaining at least one experimental result from an experimental scan of at 
least one Region of Interest; and 

determining at least one resonance parameter for at least one experimental 
result relative to at least one calibration profile. 

2. The method of claim 1, wherein one resonance parameter is an angle shift. 

3. The method of claim 1, wherein one resonance parameter is a wavelength 

shift. 

4. The method of claim 1, further including the step of storing at least one 
caUbration profile in memory. 

5. The method of claim 1, fiirther including the step of storing at least one 
resonance parameter in memory. 

6. The method of claim 1, further comprising the step of computing at least one 
calibration set statistic. 

7. The method of claim 6, fiirther including the step of displaying at least one 
caUbration set statistic. 
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8. The method of claim 1, wherein said step of generating a calibration profile for 
at least one scanned Region of Interest comprises the steps of: 

generating a raw calibration profile; and 

determining at least one derivative of said calibration profile from the raw 
calibration profile. 

9. The method of claim 8, wherein said step of generating a calibration profile for 
at least one scanned Region of Interest further con5)rises the step of smoothing said raw 
calibration profile. 

10. The method of claim 8, wherein said step of generating a calibration profile 
for at least one scanned Region of Interest further comprises the step of determining at 
least one property of said calibration profile from the raw calibration profile. 

1 1 . The method of claim 10, wherein the properties determined are selected from 
the group consisting of Full Width at Half Maximum, nominal resonance angle, fractional 
depth, and maximum intensity. 

12. The method of claim 9, wherein said step of generating a calibration profile 
fiirther comprises the step of subsampling the smoothed raw calibration profile. 

13. The method of claim 12, wherein said step of generating a calibration profile 
further comprises the step of extrapolating the ends of the subsampled smoothed raw 
calibration profile. 

14. The method of claim 1 3, wherein said step of generating a calibration profile 
further comprises the step of performing a second smooth of the subsampled smoothed 
raw calibration profile. 
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15. The method of claim 14, wherein said step of generating a calibration profile 
further comprises the step of storing the calibration profile in memory. 

1 6. The method of claim 15, wherein said step of generating a calibration profile 
further comprises the steps of: 

determining the quality of the calibration profile; and 

marking the calibration profile according to the quality determination. 

17. The method of claim 1, further including the step of performing a preliminary 
quality check on at least one calibration result. 

18. The method of claim 17, further including the step of flagging at least one 
calibration result in memory as valid or invahd according to the results of the preliminary 
quality check. 

19. The method of claim 1, further including the step of computing a derivative 
of at least one calibration profile. 

20. The method of claim 1, further including the step of displaying at least one 
scan result to a user. 

21. The method of claim 1, wherein said step of determining at least one 
resonance parameter for said experimental scan of at least one Region of Interest 
comprises the steps of: 

calculating an estimated resonance shift; 

calculating at least one interpolated profile firom said estimated resonance 
shift and said calibration profile; 

fitting said experimental scan, using said interpolated calibration profile; 

obtaining fit coefficients from said step of fitting; 

calculating, fit)m the fit coefficients, a residual resonance shift from the 
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resonance shift; 

calculating an improved estimate of the resonance shift; and 
iterating until the value of the resonance shift converges to a 
predetermined convergence criterion. 

22. The method of claim 21, wherein said step of determining at least one 
resonance parameter for said experimental scan of at least one Region of Interest 
comprises the step of calculating fit residuals. 

23. The method of claim 21, wherein said step of determining at least one 
resonance parameter further includes the step of estimating the time of scan minimum. 

24. The method of claim 23, wherein said step of determining at least one 
resonance parameter further includes the step of initially pruning the experimental scan to 
within the Umits of the calibration profile. 

25. The method of claim 24, wherein said step of determining at least one 
resonance parameter further includes the step of fitting to a sweet zone, comprising the 
steps of: 

truncatrug the interpolated profile to the sweet zone; and 

redetermining the resonance parameter utilizing the truncated interpolated 

profile. 

26. The method of claim 25, wherein said step of determining at least one 
resonance parameter further includes the step of performing initial data validity checks. 

27. The method of claim 26, wherein said step of performing initial data validity 
checks comprises the steps of: 

checking profile availability; 
checking self-consistency of data; and 



-59- 



• ' PATENT 



checking scan indexing. 

28. The method of claim 21, wherein said step of fitting employs a least squares 

fit. 

29. The method of claim 1, further including the step of reporting errors in an 
error log. 

30. The method of claim 29, wherein the step of reporting errors utilizes a local 
error log. 

3 1 . The method of claim 29, wherein the step of reporting errors employs remote 
error reporting. 

32. The method of claim 1, further including the step of performing a chip 
qualification check. 

33. A method for quantitation of surface-binding optical resonance profiles 
comprising, in combination, the steps of: 

obtaining at least one calibration result firom a calibration scan of at least 
one Region of Interest; 

generating, firom at least one calibration result, a caUbration profile for at 
least one scanned Region of Interest, comprising the steps of: 
generating a raw calibration profile; 
smoothing said raw calibration profile; 
subsampling the smoothed raw calibration profile; and 
determining properties of said caUbration profile firom the 
smoothed raw calibration profile; 
storing at least one calibration profile in memory; 
computing a derivative of at least one calibration profile; 
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obtaining at least one experimental result from an experimental scan of at 
least one Region of Interest; 

determining a resonance shift of at least one experimental result relative to 
at least one calibration profile, comprising the steps of: 
calculating an estimated resonance shift; 
calculating at least one interpolated profile from said estimated 
resonance shift and said calibration profile; 

fitting said experimental scan, using said interpolated calibration 

profile; 

obtaining fit coefficients from said step of fitting; 
calculating, from the fit coefficients, a residual resonance shift 
from the resonance shift; 

calculating an improved estimate of the resonance shift; 
calculating fit residuals; 

iterating until the estimated value of the resonance shift converges 
to a predetermined convergence criterion; and 
estimating the time of scan minimum; and 
displaying at least one scan result to a user. 

34. An apparatus for quantitation of surface-binding optical resonance 
profiles comprising, in combination:.. 

calibration module, said calibration module conq)rising: 

calibration scan result fetcher; and 

calibration profile creation module; and 
fitting module, said fitting module comprising: 

experimental scan result fetcher; 

calibration profile fetcher; and 

resonance shift determination module. 

35. The apparatus of claim 34, wherein said calibration profile creation module 
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further includes a curve smoother. 

36. The apparatus of claim 34, wherein said calibration profile creation module 
further includes a subsampler, 

37. The apparatus of claim 34, wherein said calibration profile creation module 
further includes a curve smoother and a subsampler. 

38. The apparatus of claim 34, further including a resonance parameter 
calculator. 

39. The apparatus of claim 38, wherein the calculated resonance parameters are 
selected from the group consisting of estimated absolute resonance point, time of 
resonance minimum, diagnostic information, and quality information. 

40. The apparatus of claim 34, further including an instrument control and data 
acquisition module. 

41 . The apparatus of claim 40, further including a test and support module. 

42. A method for quahfying a surface plasmon resonance chip comprising, in 
combination, the steps of: 

obtaining a golden calibration profile for the type of chip to be qualified; 

obtaining at least one calibration result from a calibration scan of at least 
one Region of Literest of a chip to be tested; 

comparing said at least one calibration result to said golden calibration 
profile to obtain at least one comparison result; and 

determining whether said chip is suitable for use by applying selection 
criteria to said at least one comparison result. 
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43. The mettiod of claim 42, further including the step of displaying chip 
qualification results to the user. 

44. The method of claim 43, wherein the step of determining whether the 
chip is suitable includes the step of incrementing a "bad ROI" count. 

45. The method of claim 44, wherein the step of determining whether the 
chip is suitable includes the step of storing a "bad ROI" number for display. 

46. The method of claim 42, wherein flie step of comparing the 
calibration includes the step of initializing a fit module with a chip qualification 
parameter set. 
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ABSTRACT 

Empirical Profile Fits are used to quantitate Surface-binding Optical Resonance 
profiles. The EPF process has two stages, Calibration and Fit, In the Calibration stage, a 
calibration surface-binding optical resonance scan is obtained with relatively fine angle 
5 or wavelength spacing over a range sufficient to include fiill resonance profiles for all 
regions. Smoothed, subsampled empirical profiles for each Region of Interest are 
generated by the main caUbration module, together with first derivative curves and 
diagnostic information. Properties returned may include approximate resonance position, 
depth, and width. In the Fit stage, individual ROI scans are used for measurement of 

10 resonance shifts relative to the empirical profile. The fitting module identifies the region 
of the experimental scan encompassing the resonance and fits that region using the 
previously stored empirical profiles, quantifying and returning the desired values, 
including the shift in the resonance as compared to its location at calibration, estimated 
absolute angles or wavelengths, time of resonance minimum, and additional diagnostic 

15 and quahty information. Optionally, data obtained fi-om either Calibration or Fit stages 
may be exported for analysis on oflier systems. In a preferred embodiment, Instrument 
Control and Data Acquisition Software sets intemal parameters in the EPF Calibration 
module and sends raw data firom a calibration scan to the EPF Calibration module, which 
funnels the data through a Subsampler and a Savitsky-Golay smoothing routine before 

20 taking derivatives and characterizing the data to create an empirical profile for the chip. 
The empirical profiles are then optionally stored. Next, Instrument Control and Data 
Acquisition Software sets intemal parameters in the EPF Fitting module and sends raw 
data firom a run-time scan performed utilizing the chip to the EPF Fitting module, which 
qualifies it, queries the EPF Calibration module for the empirical profile for the chip, and 

25 fits the curve, iterating when necessary. Results firom the fitting process are then returned 
and thereby provided to the user. 
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